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ABSTRACT 


This  report  presents  results  of  a  study  on  elastic-airframe  dynamics  that 
are  in^jortant  froir<  the  standpoint  of  flight  control  system  design.  Approximate 
transfer  functions  are  given  In  literal  terms  for  three  classes  of  vehic3.es. 
These  ai-e  of  such  a  form  that  the  important  poles  and  zeros  are  related  directlj’ 
to  simple  function-  of  aerodynamic,  elastic,  and  inertial  properties.  The  aero- 
elastic  corrections  required  to  account  for  the  flexibility  influences  of  all 
modes  not  included  in  the  equations  of  motion  are  discussed,  and  a  rigorous 
method  for  applying  these  corrections  is  presented. 
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SLCTION  I 


IMTRODUCTION 


A.  GENERAL 

The  high-speed  capabilities  of  modem  airplanes  depend  on  (among  other  things) 
the  use  of  extremely  low  thickness  ratios  for  lifting  surfaces  and  on  very  high 
fineness  ratios  for  bodies.  Coupled  with  desired  payload  and  range  capabilities, 
which  impose  natural  restrictions  on  weight,  this  dependence  leads  to  fairly 
flexible  structures  and  relatively  low  frequencies  for  the  structural  oscillatory 
modes.  For  certain  flight  conditions,  these  modes  tend  to  couple  with  the  rigid- 
body,  short-period  motions;  in  some  cases,  this  tendency  is  greatly  exaggerated 
by  the  action  of  the  autopilot.  The  danger  of  such  autopilot-flexible  airframe 
coupling  generally  increases  as  the  structure  is  lightened  to  reflect  reduced 
stiffness  requirements.  In  such  cases,  the  incipient  airframe-autopilot  instabil¬ 
ity  must  be  checked  by  analyses  which  may  require,  in  addition  to  the  normal 
rigid-body  modes,  consideration  of 

1..  as  many  as  the  first  three  or  four  coupled  normal  ("free-free") 
modes,  which,  in  general,  comprise  fuselage,  wing,  and  tail 
deflections 

2.  structural  damping  effects,  usually  included  as  an  equivalent 
viscous  damping 

3.  contributions  of  the  structural  modes  to  the  sensor  output 

4.  control  system  nonlinearities  and  more  detailed  treatment  of 
control  system  dynamics  than  necessary  for  rigiu  situations 

5.  nonstationary  aerodynamic  effects. 

Regardless  of  its  complexity,  the  closed- loop  system  must  be  stable  for  each 
of  the  flight  regimes  to  be  encountered.  Additionally,  it  must  accept  the 
required  guidance  inu"*:,  .-...d  'ual  with  ur,.lc«-'lred  inputs,  such  as  atmospheric 
turbulence  and  noise  generated  by  the  airframe-autopilot  sys+em  itself,  and  must 
not  exceed  structural  or  other  limits-  Preferably,  this  is  to  be  accomplished 
with  a  6in5)le  control  system. 


Manuscript  released  by  the  authors  April  1962  for  publication  as  an  .ASD 
'i’ci'.nnic?!  I'ocumer.tar:'  P.oporl 


ASD-TDR-62-279 


1 


M  hoc  solutions  to  the  problem  posed  above  liave  been  obtained  by  increasingly 
complex  multidegree-of ‘freedom  analyses  which  involve  the  use  of  large  computer 
facilities.  Such  analyses  provide  little  in.  .jiit  into  the  physics  of  the  problem j 
consequently,  interpretation  of  the  results  to  obtain  more  than  a  yes-no  answer 
and  an  extension  ''f  the  findings  to  slightly  modified  situations  is  difficult. 
Furthermore,  there  is  little  carry-over  from  system  to  system,  so  that,  for 
example,  the  mjmber  and  types  of  basic  degrees  of  freedom  required  to  yield  the 
critical  situation  for  a  new  design  cannot  readily  be  assessed  a  priori.'  The 
design  process  suffers  accordingly.  Not  only  is  an  undesirably  long  time  requii-ed 
to  select  (and  perhaps  later  modify)  the  pertinent  degrees  of  freedom,  and  to  set 
up  and  run  the  problem,  but  also  modifications  required  to  cure  discovered  problem 
areas  are  difficult  to  explain  to  a  design  group  affected  or  to  management.  For 
these  reasons,  it  is  desirable  to  obtain  simple  literal  approximations  to  the 
airframe  transfer  function  factors.  Such  approximations  relating  the  important 
poles  and  zeros  directly  to  simple  functions  of  aerodynamic,  elastic,  and  inertial 
properties  can  provide  an  invaluable  design  guide. 

This  report  presents  the  results  of  a  1-year  stu<3y  devoted  to  the  analytical 
approximation  of  flexible  airframe  transfer  functions.  In  this  study,  three 
classes  of  vehicles  were  represented  by  typical  configurations,  and  the  influence 
of  elastic  modes  on  the  longitudinal  transfer  functions  was  examined  (the  three 
configurations  were  also  subjected  to  a  parallel  study,  Ref.  i,  which  yielded 
the  basic  Information,  i.e.,  mode  shapes,  etc.,  used  in  this  study).  Each  of 
the  cord’igurations  is  shown  to  possess  transfer  function  facfors  which  can  be 
simply  approximated  by  a  few  major  terms  when  interest  is  confined  to  only  the 
first  two  elastic  modes.  (The  x'esults  of  Ref.  1  show  that  additional  higher- 
frequency  modes  do  not  appreciably  affect  the  two-elastic-mode  transfer  fimctlon 
in  the  frequency  range  that  is  important  for  flight  control  analysis.) 

The  first  configuration  studied  (Configuraticni  2)  is  a  missile-like  vehicle 
with  canard  control,  capable  of  supersonic  speeds  at  low  altitudes.  Configura¬ 
tion  3  is  a  swept-wing,  high-aspect- ratio  arrangement,  while  Configuiation  4  is 
a  supersonic  delta-wing  vehicle.  All  three  configurations  are  described  in 
detail  in  Appeniix  C. 
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B,  OUTLINli  OF  THE  REPORT 


The  report  i-.s  divided  into  five  Bcctions  ^  Mowed  by  four  appendices.  Most 
of  the  analytical  work  is  presented  in  Sections  II  through  IV,  while  the  numeri¬ 
cal  data  are  included  in  the  appendices. 

Section  II  presents  a  derlration  of  the  equations  which  are  used  to  form  the 
transfer  functions.  The  static  aeroelastic  effect  of  the  truncated  high-l'vequency 
modes  on  the  aerodynamic  inputs  Is  discussed  and  the  method  of  inclusion  is  shown. 
The  derivation  of  the  matrix  required  in  this  method  is  presented  in  Appendix  A. 

The  transfer  function  factored  forms  are  presented  in  Section  III,  along  with 
the  approximation  formulas  for  the  three  configurations  studied.  All  of  these 
approximations  were  derived  by  one  of  the  methods  presented  in  Appendix  B. 

Section  IV  discusses  the  problem  of  sensor  location  for  closed-loop  operation. 
The  effect  of  sensor  location  on  a  particular  configuratioi  is  shown,  and  a  method 
for  "optimum"  placement  is  suggested.  ("Optimum"  here  Implies  that  the  effect  of 
the  elastic  modes  is  minimized  with  respect  to  the  rigid-body  pitch  degree  of 
freedom.) 

As  noted,  ^pendix  A  presents  a  derivation  of  a  proper  method  of  accounting 
for  the  elastic  modes  not  included  in  the  equations.  This  method  was  derived  in 
Ref.  1 ,  and  Appendix  A  parallels  that  presentation.  An  exaiiq>le  is  included  which 
utilizes  the  method,  and  shows  the  exactness  of  the  results  obtained. 

Appendix  B  presents  the  several  methods  that  were  used  to  derive  the  transfer 
function  approximation  fonnulas  of  Section  III.  No  one  method  could  be  found 
which  consistently  produced  the  simplest  approximations;  hence,  the  approximations 
were  derived  by  the  best  of  those  in  Appendix  B  for  the  case  at  hand. 

A  detaided  description  of  the  configurations  studied  is  presented  in  Appendix  C, 
along  with  the  noiinal  mode  shapes  used. 


The  numerical  equations  of  motion  for  each  conl’lguration  are  included  in 
Appendix  D,  as  are  the  exact  transfer  functions  that  these  equations  yielded. 
The  equations  and  the  transfer  functions  were  calculated  by  a  digital  computer 
according  to  the  equations  outlined  in  Section  II.  The  transfer  functions 
otJi....-r'''d  with  the  approximatsoi.  foimulas  of  Section  III  are  also  presented  in 
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Appendix  D  as  an  indication  of  their  accuracy.  All  of  these  data  are  presented 
for  a  range  of  dynamic  pressures  for  each  configuration  with  one  and  two  elastic 
inodes  included. 
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SECTION  II 


EQUATIONS  OF  MOTION 

The  eq,aations  will  first  be  derived  in  terms  of  the  physical  coordinates  of 
the  airframe,  and  will  then  be  converted  to  modal  coordinates  to  allow  a  reduc¬ 
tion  in  degrees  of  freedom  and  coupling  tenns  by  the  use  of  orthogonal  modes. 

The  high-frequency  modes  will  then  be  eliminated  and  the  equations  reduced  to  a 
set  involving  a  limited  number  of  flexible  modes.  Detailed  consideration  of  the 
anvodynamic  forces  will  then  give  the  form  of  various  coefficients  involved  in 
che  final  equations  of  motion. 

A.  GENERAL  EQUATIONS  OF  MOTION 

The  methods  for  the  development  of  eypressions  for  the  inei’tial  and  elastic 
forces  on  a  flexible  airframe  differ,  depending  on  whether  the  inertial  and 
external  loadings  are  considered  to  be  distributed  or  concentrated.  If  they  are 
considered  to  be  concentrated  at  a  finite  number  of  points,  then  the  displacement 
of  any  point,  qj^,  can  be  written 

Qi  =  E  CijQj  (l) 

where  q^^  ir  the  physical  displacement  (which  may  be  either  a  translation 
or  n  rotation)  at  the  i^-h  point 

C^j  is  the  flexibility  influence  cor?fficient  givxi.ij  the  phj'Cicai 
displacement  at  the  ith  point  caused  by  a  unit  physical  force 
(i.e.,  a  force  or  a  moment)  at  the  jth  point 

Qj  is  the  physical  force  applied  at  the  jth  point 

A  force  and  a  moment  o"n  >■?  simultaneously  at  any  given  location  merely 

by  making  two  of  the  n  points  of  application  coincident;  e.g.,  th'  first  and 
second  of  the  n  points  will  be  the  same  if  Q;  is  a  force  applied  at  some  locatic-, 
and  Qg  is  a  moment  applied  at  the  same  location. 

The  physical  elastic  force  at  point  i  due  to  any  arbitrary  set  of  ph.Y8lcal 
displacements,  qj,  can  be  expressed  as 
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n 

Qci  =  -  p  kijqj 

J=1 


(«) 


where  is  a  stiffness  influence  coefficient.  For  any  j  (e.g.,  j  =  r),  the 
set  of  is  equal  to  the  set  of  forces,  Q^,  required  to  make  q^.  =  1  and 

qj^  =  0  for  i  r. 


If  the  distributed  air  loads  over  the  airframe  are  considered  to  act  as  a  set 
of  concentrated  forces,  the  sum  of  the  elastic  and  aerodynamic  forces  at  the  i^h 
niass  may  be  equated  to  the  inertial  force  at  that  point,  yielding  (neglecting 
structural  damping) 


mi 


«i^qi 

dt2 


Qe 


.  +  a,. 

i  “1 


(3) 


where  m^  is  the  physical  mass  at  the  i^h  point 

is  the  physical  aerodynamic  force  at  the  i'^h  point 

The  aerodynamic  force  may  be  considered  to  be  generated  by  the  displacement  of  a 
finite  number  of  physical  coordinates  on  the  airframe. 


Qai  =  E  Rir*J 

J  — I 


(4) 


where  is  the  physical  aerodynamic  force  at  point  i  caused 

by  a  unit  movement  of  point  J.  R^j  will  generally 
be  a  polynomial  in  the  differential  operator  (or 
Laplace  variable),  s 


Using  matrix  notation,  it.  Is  nossib.le  to  write  2q  2  and  4  as 


{%}  = 

-Mlat 

(5) 

M  kl 

(6) 
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Equation  3  now  becomes 


MW  .  -Mk’  MW  (7) 

It  '.s  customary  to  combine  the  elastic  and  inertial  forces  because  these  do  not 
vary  with  dynamic  pressure: 

[ras2+k]{q|  =  [R]{q}  (8) 

If  movement  of  the  control  surfaces,  while  introducing  aerodynamic  forces  into 
the  system,  does  not  introduce  significant  inertial  or  elastic  force,  it  is 
reasonable  to  separate  control  surface  deflections  from  the  rigid-body  and  elastic 
deflections,  and  to  write  Bq  8  as 

[ras2  +  k]{qj  =  [R]{q|  +  {Qin}  (9) 

where  is  the  physical  aerodynamic  force  input  caused  by  the  move¬ 

ment  of  control  surfaces.  Gust  loads,  nonuniform  wind 
conditions,  etc.,  will  create  forces  which  can  also  be 
included  in  the  term 

Equation  9,  when  expanded,  appears  as  follows; 

mis2q,  +  k^^q,  +  k^gqG  +  +  k,j,qn  =  R,.,q,  +  R,2q2  +  •"  +  R,n9n  '^ini 

mos^qg  +  k2i  q^  +  k22q2  +  •  •  •  +  k2nqn  =  +  ^22^  +  •  •  •  +  R2nqn  ^n2 

(10) 

•  • 

«  • 

“ns^%  +  +  K2^2  +  •••  +  kn^qn  =  R^.q,  +  Rngqg  +  *••  + 

Equations  10,  then,  are  the  equations  of  motion  (ioi  perturbations  from  a 
trimmed  condition)  for  a  flexible  airframe  represented  by  n  control  points  (see 
Fig.  1).  The  degree  of  accuracy  employed  in  the  construction  of  the  [r]  matrix, 
and  the  number  and  location  of  discrete  mass  points  chosen,  will  determine  the 
adequacy  of  these  equations  in  representing  the  actual  system.  Because  a  great 
nuE'her  (  f  q, ’ s  io  g-Hioral-ly  ,..J  red  for  an  adequate  representation  of  a 
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continuous  airframe  with  a  contlnuoun  loading,  It  Is  lii5)ractlcal  to  work  directly 
with  Eq  10  to  achieve  simple  approximate  methods.  As  an  equation-reducing 
alternative,  airframe  motion  is  often  represented  by  a  few  normal  modes;  then, 
each  of  the  qj^'s  consists  of  the  superposition  of  the  motion  nf  the  modal  coordi¬ 
nates;  l.e., 

qi  =  g  (11) 

J“  I 


Defining 


(12) 


Then, 

n 

Defining 

Then, 

n 

11 

(13) 

(111) 

(15) 


where  q^  is  the  physical  displacement  at  the  il-h  point  on  the  airframe 

(ft.  is  the  physical  displacement  of  the  i'tli  point  caused  by  a 
“  unit  generalized  displacement  of  the  jth  normal  mode 

is  the  physical  normalized  displacement  at  the  i'th  point 
caused  by  a  unit  generalized  deflection  of  the  jl'h  normal 
mode.  The  collection  of  all  the  (pj^'s  for  any  given  j 
reprsseiii-B  the  shape  for  the  jth  normal  mode. 
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J  j  is  the  generalized  displaoerneat  or  coordinate  of  the 

normal  mode;  i.e.,  is  a  scale  factor  for  the  normal 
mode,  g^ven  by  the  normalized  physical  displacement  (result¬ 
ing  from  deflection  of  the  normal  mode,  and  no  other)  of 
a  preselected  point  on  the  ai”^  ■ne. 

q,  =  +  9,262  +  •••  +  <P,n6n 

52  =  ‘l>2ih  ^Sh  +  •••  +  ‘P2n5n 

h6) 

* 

5n  =  ‘Pnl^l  +  <Pn252  +  •*'  +  'PnnSn 


This  may  be  written  in  matrix  form  as 


lq|  -  M)S( 

(17) 

Henceforth,  [$]  will  be  referred  to  as  the  modal  matrix, 
presentation  of  Eq  16  is  given  in  Fig.  2. 

A  typical  graphical 

The  n  columns  of  [<1>],  l.e.,  the  mode  shapes,  are  found  by  assuming  siinple 

harmonic  motion  (s  =  jo)  and  substituting  Eq  17  into  Eq  9, 

[k  -  mu|]{9i}  =  {0} 

(18) 

which  can  also  be  written 

Mi<Pi}  -  [®4ll<Pi}  =  0 

(19) 

[kj{9i^  -  =  0 

(20) 

tm]"’[k]{9j  =  4|9ij 

(21) 

The  form  of  Eq  21  makes  it  clear  that  the  o^'s  are  eigenvalues  of  ,  and 

the  9j^'s  are  the  associated  eigenvectors. 

Foi‘  any  two  linearly  independent  eigenvectors,  it  is  possible  to  write 
(Eq  20) 
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n 

Firwi’1,2.  riocl';  CiorvrG;  <ij^  =  21 

j-1 
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(22) 


DO{<J>a}  =  “In  {'Pal 

MKl  •=  “|NW  (25) 

Vfhen  both  sides  of  Eq,  22  are  transposed, 

a'hereforsj 

But  [k]  and  [n^  are  symmetric,  whereby 

kaJW{<Pb}  =■  “IbflJ  N{<Pbl  (24) 

Now,  premultiplying  both  sides  of  Eq  25  by  , 

L<pUMK}  =  “|L<€jNKf  (25) 

A  comparison  of  Eq  24  and  25  yields 

“I  N  I'Pb}  =  “I  LfU  N  i'Pbl 

or 

(“I  -  “|)  L<p|l  N{<Pb}  =  0  (26) 

If  a:|  and  oj^  are  distinct,  then 

b|lN{<Pb}  =  0  (27) 


Section  1 .21  of  Ref.  5 
shows  that  Eq  2J  above 
off-diagonal  terms  are 


considers  the  case  where  and  are  "'^t  distinct,  and 
still  holds.  Since  the  implication  of  Eq  27  is  that  all 
zero,  then  for  orthogonal  modes 


where  must  be  diagonal.  It  can  similarly  be  shown  that 


.  M  (2?) 

where  must  be  diagonal. 

The  significance  of  Eq  28  and  29  is  that  Eq  9,  combined  with  Eq  17,  can  be 
premuitiplifed  by  [''j  to  obtain 

=  [■»]*[»]  H(!|  +  (30) 

where  [4']*[ms^  +  k]  [<i>]  will  be  diagonal.  Therefore,  defining 

[Ms2  +  k]  h  [<t.]^[ms2  +  k]  [c] 

[P|]  =  bfim  (5-) 

|Fin^  = 

allows  Eq  JO  to  take  the  simple  form 

[ms2  +  k]{|{  =  EF|]il}  +  {^inl  (52) 

represents  the  generalized  mass  matrix  in  generalized 
coordinates 

represents  the  generalized  stiffness  matrix  in  generalized 
coordinates 

is  a  column  of  generalized  coordinates  (which  ai-e  orthogonal 
coordinates  when  [F|3  =  O) 

represents,  in  generalized  coordinates,  the  externally 
applied  forces  per  unit  deflections  in  |||  (and  unit 
velocities  in  {1}  ,  etc.) 

The  left  side  of  Eq  J2  represents  the  structural  dy.'iamics  of  the  vehicle  in 
vacuum  (neglecting  structural  damping),  while  the  right  side  represents  externally 
applied  forces.  When  expanded,  the  left  side  will  appear  as 


where  [mJ 
[K] 
{i} 
[F§] 


(g^  + 

M2(s2  +  u|)62 

<  .  > 


As  is  normal,  I-;  and  5 2  will  be  used  to  represent  rigid-body  translation  and 
rotation,  respectively.  Thus, 

=  m,  the  total  physical  mass  of  the  system 
M2  =  ly,  the  total  physical  pitch  inertia  of  the  system 

=  h  (positive  down)  (55) 

I2  =  6  (positive  nose  up) 

<0|  =  =  0  (no  structural  stiffness  in  li  or  62) 


The  appearance  of  h  is  not  in  keeping  with  the  normal  aircraft  stability  and 
control  formulation  of  the  equations  of  motion.  Therefore,  it  is  desirable  to 
transform  Eq  52  so  that  h  is  replaced  by  w.  Such  a  transformation  is  represent¬ 
able  by  the  matrix  [w] ,  where  jw]  is  defined  by 


''h " 

s 

e 

e 

\  ^ 

; 

/ 

♦ 

• 

• 

V  •  > 

• 

1  •  y 

(5*^) 


1  V  iiiih  case,  the  matrix  [wj  can  be  found  rather  easily  by  using  the  expression 
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relating  v  to  h  and  6  for  the  assumed  linperturbed  condition,  =  0. 


Figure  5*  Derivation  of  the  (w]  Matrix 


Thus, 


and  the  desired  transformation  matrix  is 


1 


s 


0  0 


0 


0  10  0 


0 


[w] 


0 

0 

0 


0  10  0 

0  0  1  0 

•  *  % 

•  •  • 

0  0  0  ...  1 


(55) 


Using  the  transl'orm.'.i-dcc  mat.ix  deij.;.^.d  cy  Sti  5^,  rnd  for  convenience  defining 


[y]  =  [ms2  + 

Eq,  52  becomes 

MMiit  =  WMW  ♦  Knt 


(56) 
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where  ||  is  now  w/s  rather  than  h.  It  can  be  shown  that  [wj  is  not  a  diagonal 
matrix.  It  is  therefore  desirable  to  redefine  to  include  the  off-diagonal 

term  from  MM,  thus  giving  a  diagonal  form  to  M  [w] .  This  can  be  done  by 
writing  out  the  MCw]  matrix, 


M  M 


ms^  -UQms  0 

0  lyS^  0 

0  0  MjCs^  +  a^) 


0 

0 

0 


(38) 


0 


0 


0 


kf  /-r 

‘*nv" 


€md  redefining  [F|]  [W]  so  that  it  includes  the  -Uoms  term  that  appears  in  E(i  58. 
This  will  restore  M  [w]  to  its  original  diagonal  form. 

B.  ELIMINATION  OF  THE  HIGH-FREQUENCY  MODES 

At  this  point  the  set  of  simultaneous  relationships  given  by  Eq  37  Is  capable 
of  yielding  results  which  increase  in  "exactness"  with  the  number  of  modes  con¬ 
sidered.  Since  engineering  interest  is  inevitably  confined  to  a  limited  band¬ 
width,  the  importance  of  including  higher  frequency  modes  ic  measured  by  their 
effects  in  this  bandwidth.  For  example,  if  all  elastic  modes  are  considered  to 
lie  outside  the  frequency  region  of  concern  and  all  are  excluded  from  the  equa¬ 
tions,  then  the  resulting  solution  yields  only  the  conventional  rigid-body 
short-period  motions.  But  at  appreciable  dynamic  pressures  this  is  a  gross  over¬ 
simplification,  becau.-c  chc  ex^ludea  1..  .laclural  give  rise  then  to  at  least 

an  aeroelastic  correction  on  the  rigidrbody  stability  aeiivalivcs.  bucn  correc¬ 
tions  can  be  and  usually  are  made  by  considering  only  the  static  deflection 
properties  of  the  structure.  In  this  instance,  the  effects  of  all  possible 
elastic  modes  have  been  approximated,  for  the  frequency  region  of  interest,  by 
considering  only  static  deflection  characteristics.  When  the  bandwidth  of  interest 
in-.lude-  H  number  of  iov.--fr.''j>'.‘.jcy  structural  modes,  the  question  arises  as  to 


the  proper  aeroelaotio  correction  whereby  to  approxlmte  the  influence  of  ths 
neglected  modes.  Clearly,  now  the  use  of  the  "full"  aeroelastic  correction  will 
be  incorrect  since  the  modes  included  in  the  equations  of  motion  must  somehow 
alter  the  approximated  contribution  of  all  remaining  modes.  The  proper  treatment 
of  the  neglected  higher-f.requency  modes  first  studied  in  Ref.  1  will  now  be 
outlined. 

If  Eq  32  is  partitioned  into  those  coordinates  which  are  of  interest  (denoted 
5o+f>  where  o+f  stands  for  zero  frequency  +  finite  frequency),  and  into  those  of 
higher  frequency  which  are  not  of  direct  interest  (denoted  as  Iqq ) >  then  the 
latter  can  be  eliminated  from  the  equations.  This  is  accomplished  as  follows: 
Using  Eq  3I  and  3^,  Eq  32  is  rewritten  as 

Mlh  *  W^MWUI  +  W'toint  (») 

where  the  transformation  from  =  h  to  S-j  =  w/s  has  not  yet  been  made.  Then, 
partitioning  the  joatrices. 


and  expanding  the  right  side 


('>1) 


. . .  i  * '  ‘  *  I  * ' '  ’  ^ 
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Equation  is  equivalent  to  the  two  simultaneous  matrix  equations. 


tXo+^Uo+f}  =  +  CWl^MC'I’colUrt?}  +  C'^oifl’^{%n}  (^*2) 

and 

MhJ  =  [♦»fM[w]{w}  +  [U’'Me*Jt6oo)t&cJ®W  W) 

Equation  43  can  be  solved  for  ISqq},  which  can  then  be  substituted  into  Eq  42. 

This  will  result  in  equations  of  motion  for  the  modal  coordinates  of  interest. 
Thus,  multiplying  Eq  43  by  [YqJ  and  solving  for  {^oo} 

i^ool  =  [M  “  &oo]  [iIIqq]  [r]  [^oo]^  [Yqo]]  ^  OfooD^Dil  [‘‘D+fllSo+f} 

*  [m  -  W'[»oo]’'[ia:w]'’&»r'[»»:^kJ 

(4lt) 

and  using  this  result  in  Eq  42, 


[Yo+f]{lo+f}  =  Cwl'^WCWllSo+f} 

+  C«'o+f]’^[R]  C-fcJ  [M  -  Otool]  'b'ad"’B>a>f[R]&o+d 

-  CWD'^CR]  W[W  -  I^]t*cx]"[R]&co]I’'^"&a>f|%„} 

"  (43) 


Because 


U[d][b]  +  [d][c]]  ’  =  [[b]  +  [cj]’^[d]-’ 
[MCbJ  +Cc][b]]''  =  [b]-'[[d]  +[c]]'' 
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Then 


[W  - 

»  -  W  [•>„]]■' 

Using  this  result,  Eq  45  is  modified  to 

C^o+d  { 6o+f }  =  C'*>o+f]  ^ to  O^o+f]  Uo+f } 

*  -  MWM-’&jf  WCwlliow’) 

» C*««f61WiyM|iil  -  MWW wf’ Knt 

*  fet/ldinl 

and  collecting  terms, 

Dfc+fliSo+f]  =  Cwl'^ 

«  [m  *  MM &«]■'[♦  j’'[m  -  mmm-'mT] 

X  [[R][w3Uo+ff  +  {Qin}]  (^7) 

To  simplify  this  further,  use  the  following  identity: 

[l]  +  [b]  [UJ  -  [b]]  '  .  [I]  ♦  [b]  [[i]  -  [B]]  '  .  [M  -  [B]]  '  +  [[I]  -  [3]]  ' 

-  W  -  [W  -  Lb]][li]  ■  tB]]’’  ^  [[I]  -  [Bj]'' 

-  H  -  M  ♦  [W  -  [»]]■' 

-  [M  -  Cb]]'' 
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whereby  Eq  47  reduces  to 


«o.,  dtflne  [A]  .  [h  -  MM  (W 

Equation  48  then  becomes 

[Yo+^IwJ  =  CW]''rA][R][W]IW}+  CWl^CA]Kn}  (50) 

Equation  50  Is  seen  to  be  very  similar  to  Eq  59-  However,  In  Eq  50,  the 
coiJmns  of  modal  coordinates  contain  only  those  coordinates  which  are  of  Interest j 
therefore,  this  equation  represents  a  fewer  number  of  sim) Itaneous  differential 
equations  to  be  solved,  (in  essence,  the  last  r  equations  have  been  used  to 
eliminate  the  last  r  varinbles  from  the  set  of  n  simultaneou3  equations.)  Also, 
a  new  term  appears  in  Eq  50  which  was  not  found  in  Eq  59,  this  term  is  the  [a] 
matrix.  It  represents  the  modifications  which  must  be  made  in  the  first  n-r 
equations  to  include  the  effects  of  the  higher- frequency  modes.  The  [a]  matrix 
thus  represents  an  aeroelastic  correction  factor  to  the  system.  It  will  theoreti¬ 
cally  account  exactly  for  all  influences  of  the  higher- freqiiency  modes.  However, 
the  exact  calculation  of  [A]  requires  all  the  information  contained  in  a  coiiq)lete 
set  of  n  equations  and  involves  an  unwieldy  inversion  of  a  mi-trix  containing 
terms  in  s  and  s^  (see  Eq  49).  At  the  frequencies  of  interest,  u),  which  are 
always  much  smaller  than  the  higher  mode  eigenvalues,  cqt,  by  definition,  the 
s  and  s2  terms,  relative  to  the  stiffness  term,  are  proportional  to  and 
(aj/<%)^,  respectively;  thus  they  can  be  neglected.  Even  then,  calculation  of 
[a]  from  Eq  49  would  :A;quirc  K-.jwleUjiA.  the  hl'U^.r-'f’requency  mode  shapes. 
Fortunately,  however,  by  neglecting  the  s  and  s^  terras,  tnereby  making  [aJ  a 
quasi-static  correction  factor,  it  is  possible  (Ref.  1 )  to  calculate  the  quantity 

■  [xj  (511 

from  a  knowledge  of  only  the  static  deflection  characteristics  of  the  Gyt.lem. 

The  dets.ils  v'’  this  calc-ilaj..u/n  are  given  in  Appendix  A. 
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C.  AERODYNAMIC  FORCES 


So  far,  no  mention  has  been  made  of  methods  whereby  the  [r]  ,  (VJ ,  and 
matrices  may  be  calculated.  Although  the  latter  two  are  not  necessarily  simple 
to  fom,  they  will  not  be  discussed  because  of  extensive  treatments  in  the 
literature  (e.g.,  R-»f.  J  and  4).  The  same  might  be  said  for  the  aerodynamic 
matrix,  [r3,  except  for  a  significant  major  difference.  While  the  [y]  and  [ij] 
matrices  may  be  derived  in  m£my  ways,  the  results  (for  a  given  physical  situation) 
will  always  be  the  same;  this  is  not  true  of  the  [r]  matrix,  which  depends 
inherently  on  the  assumptions  made  as  to  the  origin  of  aerodynamic  forces. 

For  purposes  of  the  present  study,  a  very  complete  formulation  of  [r]  is 
deemed  unnecessary,  because  it  can  only  affect  certain  of  the  numbers  appearing 
in  the  equations  of  motion,  oince  these  numbers  are  required  to  be  only'  repre¬ 
sentative  of  the  configurations  involved,  the  aerodynamic  matrix  will  be  reduced 
to  a  very  simple  form.  That  is,  almost  all  secondary  aerodyamic  and  elastic 
effects  (e.g.,  wing-body  interference,  unsteady  aerodynamics,  chordwise  bending) 
will  be  neglected;  and  the  air  forces  and  moments  will  be  represented  by  average 
derivatives  associated  with  each  lifting  surface,  or  suitable  portions  thereof. 
Accordingly,  the  center  of  pressure  (c.p.)  for  lift  is  assumed  to  be  at  a  fixed 
fraction  of  the  MAC  for  the  surface  or  portion  thereof  (0.5  for  Configurations  8 
and  4  and  0.85  for  Configuration  5):;  no  downwash  effects  are  considered  on 
Configurations  8  and  4;  the  only  (pure)  moment  is  considered  to  rooult  from 
pitching  velocity;  and  the  angle  of  attack  for  a  section  is  defined  by 


a 


(52) 


where  0  is  the  rigid  chord  rotation 

hj/ii,  is  I'.ne  -  -^leel  uisnlacement  of  the  5/4  chord 

The  selection  of  to  define  the  section  angle  of  attacK  is  in  accoinance  with 
theoretical  aerodynamics  [where,  as  a  boundary  condition,  the  flow  velocities 
over  the  upper  and  lower  surfaces  of  an  airfoil  are  matched  at  the  trailing  edge 
(Ref.  6)]. 

•Ik.  all  this  a  lit  tin  riore  specifically  to  [r]  ,  refer  to  Eq  6  where 

the  fo  are  ordered  as 


80 


.  rri 

and  the  Q's  are  ordered  as 

(Z-force  at  c.p. 
moment 


(55) 


(5'^) 


Note  that  the  moment  does  not  require  specification  of  a  point  of  application 
because  the  chord  is  assumed  rigid.  Now,  partitioning  the  [r]  matrix, 


where,  in  general. 


(55) 


a^j  is  the  Z-force  (negative  lift)  at  the  c.p.  of  the  1^^ 

"erodynamic  surface  due  to  unit  vertical  displacement  of 
th(!  3/^  chord  of  the  jth  surface 

bj^j  is  the  Z-force  (negative  lift)  at  the  c.p.  of  the 

aerodynamic  surface  due  to  unit  rigid  rotation  of  the 
chord  of  the  j^h  surface 

Cij  is  the  moment  on  the  l^h  aerodynamic  surface  due  to  unit 
vertical  displacement  of  the  3A  chord  of  the  surface 

d^j  is  the  moment  on  the  aerodynamic  surface  d’ie  to  unit 
rigid  rotation  of  the  chord  of  the  surface 


The  literal  expressions  for  the  partitions  of  [rJ  are  found  from  the  general  lift 
and  moment  equations.  The  Z-force  at  the  c.p.  (which  for  convenience  is  desig¬ 
nated  by  the  subscript  "l/4")  of  section  i  is  given  by  (excluding  downwash) 


A-,  \  . 

"T 


The  moment  on  section  i 


is  (excluding  downwash) 


(56) 


(57) 


21 


Therefore  it  is 


to  write 


matrix  are  those  comprising  [r]  and  are  given  by 

'’ij  =  (59) 

Hi  =  0 

'"ij  V 
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where  1b  the  Kronecker  b  (6  =  1  fer  1  =  j;  6^0  for  i  j),  and  the  o"her 
symbolB  are  standai'd  aerodynamic  symbols. 

The  physical  coordinates  on  the  right  side  of  Eqt  58  can  be  written 


(60) 


Therefore  Eq  5*+  becomes 


b 

d 


(61) 


and  utilizing  a  conq)atible  partitioning  of  Eq  Jl  can  be  expressed  as 


Therefore, 


(62) 


(65) 


or  in  the  terms  desired  here, 


(64) 


where 


is  the  matrix  of  mode  deflections  of  the  k-chord  points 

is  thtr  matrix  v.  mode  c-f  ■.r-.i  rigid  aerodynamic  chords, 

i.e.,  <Py  =  (d/dx)<Py,  which  is  constant  al/jng  a  ribld  chord 


Expanding  Eq  64  for  the  zero  downwash  case  yields 
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The  element  of  the  [F|]  is  I'ound  'oy  adding  the  elements  fi-om  each  of 

the  three  components. 


WC»3/du  -  ,J 

(66) 

(67) 

W[«0y  - 

(68) 

(69) 

(70) 

(71) 

-f  2  [([“tJ>,/uA5/l,J  -  (  2  “Ulfna)'  ■*  “o' 

(72) 

where  the  summation  is  over  the  n  aerodyneunic  surfaces. 

For  situations  involving  downwash,  F|^j  includes  terms  involving  off-diagonal 
elements  of  a^^j  and  (Eq  55) .  These  added  terms  are  shown  in  one  version 
the  f;.nal  equations  of  motion,  Eq  75» 

D.  FINAL  EQUATIONS  OF  MOTION 

With  all  the  important  elements  now  in  hand,  the  final  desired  equations  of 
motion  can  be  fonnulated.  To  achieve  a  form  consistent  with  aeronautical 
stability  and  control  usage  requires  application  of  the  |w|  transformation 
matrix  (Eq  5^  and  35),  litv'iwfex'i’ing  all  aerodynaraic  terms,  except  inputs,  to  the 
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left  side,  and  nondimensionalising  by  ui/iding  each  equation  by  the  appropriate 
inertial  quantity.  Note  that  the  mode  shapes  are  normalized  so  that  the  ijeneral- 
ized  mass  for  each  elastic  mode  is  unity.  Doing  this  for  the  case  where  only 
two  flexible  modes  are  included,  and  where  downwash  from  the  first  aerodynamic 
surface  affects  the  angle  of  attack  of  the  fourth  aerodynamic  surface,  results 
in  Eq  73,  whlci  utilizes  the  following  identities: 


%1 

V 


<2 


1  rigid-body  translation  imparts  equal  translation 
to  all  aerodynamic  surfaces 

0  rigid-body  translation  imparts  no  rotation  to 
aerodynamic  surfaces 

Ijj  rigid-body  rotation  imparts  translation  propor- 
^  tional  to  the  distance  from  the  center  of  rotation 
(inp.  is  the  distance  from  the  airframe  center  of 
gravity  to  the  k-chord  of  the  n^h  aerodynamic 
surface,  positive  aft) 

1  rigid-body  rotation  imparts  equal  rotation  to  all 
aerodynamic  surfaces 


It  should  be  noted  bhat  no  aeroelastic  correction,  [a]  ,  has  yet  been  applied 
to  Eq  73  and  that  to  do  so,  in  literal  terms,  would  be  a  practical  impossibility. 
Equation  73  is  thus  mainly  illustrative  of  the  form  assumed  by  the  various  coeffi¬ 
cients,  which  in  an  actual  case  would  be  modified  by  varying  aeroelastic  correc¬ 
tion  factors.  By  assigning  a  symbol  (i.e.,  a  stability  derivative)  to  each  of 
the  terms  (including  the  aeroelastic  correction),  Eq  73  can  be  simplified  and 
extended  to  the  general  situation  where  an  arbitrary  number  of  flexible  degrees 
of  freedom  are  included,  as  in  Eq  7^,  the  final  set  of  c  I'iations  of  motion. 
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SECTION  III 


FLEXIBLE  AIRFRAME  APPROXIMATE  TRANSFER  FUNCTIONS 

The  forms  for  the  longitudinal  transfer  functions  of  a  rigid  airframe  are 
well  understood,  and  a  summary  of  these  forms  may  be  found  in  Ref.  2.  The  addi¬ 
tion  of  flexible  degrees  of  freedom  to  a  system  lias  generally  been  treated  to  a 
lesser  degree,  but  the  forms  for  the  transfer  functions  are  nonetheless  also 
well  established.  In  general,  the  addition  of  each  flexible  mode  will  result 
in  the  addition  of  a  pair  of  lightly  damped  roots  to  the  numerator  and  de.nomi- 
nator  of  each  transfer  function.  Table  1  summarizes  the  forms  expected  for 
situations  where  two,  one,  or  no  elastic  degrees  of  freedom  are  included  in  the 
equations,  and  forward  speed  is  assumed  constant. 

■  In  the  current  study,  each  of  the  transfer  function  factors  shown  in  Table  I 
was  approximated  by  a  limited  number  of  terms  involving  directly  the  stability 
derivatives  appearing  in  the  equations  of  motion  (Eq  74) .  ^ese  direct  relation¬ 
ships  allow  the  effeci.s  of  parameter  changes  to  be  predicted  with  a  reasonable 
degree  of  confidence  without  actually  recalculating  the  transfer  function. 

A.  DISCUSSION  OF  METHODS  OF  DERIVATION 

Basically,  the  derivation  of  approximate  transfer  function  factors  involves 
determining  the  terms  which  are  important  for  each  airframe  configuration  con¬ 
sidered.  This  is  done  by  substituting  a  typical  set  of  numerical  values  for 
speed,  altitude,  etc.,  into  the  equations,  and  then  neglccti’".",  the  small  terms. 

In  doing  this,  it  is  assumed  implicitly  that  moderate  changes  jlii  the  parameters 
will  not  affect  the  segregation  of  small  and  large  terms;  that  is,  small  terms 
remain  small  over  a  reasonable  range  of  parameter  variation.  An  exception  to 
this  was  found  in  Configuration  4,  where  control  reversal  was  noted  for  dynamic 
pressures  of  20  psi. 

Appendix  B  contains  detailed  descriptions  of  the  two  '^i='thods  which  vf>’e  used 
to  determine  literal  approximate  factors  for  each  of  the  configurations  consid¬ 
ered  in  the  current  study.  Although  the  description  of  the  first  method  con¬ 
siders  the  case  of  factoring  a  transfer  function  denominator,  the  technique  used 
may  be  applied  to  numerators  as  wall. 
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TABLE  I 


SUMMARY  OP  TRAMSPER  FUMCTION  PACTORED  FORMS 


B.  APPROXIMATE  FACTORS 


The  approximate  factors  for  the  denominator  and  numerators  for  each  of  the 
three  configurations  are  presented  in  Tables  II  through  VII.  Inspection  of 
these  tables  reveals  that  the  transfer  function  factors  for  flexible  airframes 
contain  the  rigid  airframe  factors  d.'rived  in  Ref.  2  (with  aeroelastic  correc¬ 
tions)  along  with  tne  elastic-mode  factors.  It  is  noted  that  the  literal 
factors  for  Configurations  2  and  3  are  quite  similar,  and  in  some  cases,  are 
actually  identical.  The  results  for  Configuration  4,  however,  are  quite 
different. 

Rather  than  including  a  list  of  validity  conditions  for  each  set  of  factors, 
it  is  suggested  that  the  applicability  of  the  approximations  be  determined  by 
finding  the  exact  numerical  factors  for  a  nominal  case,  and  comparing  them  with 
the  numbers  obtained  by  using  the  approximate  formulas.  The  reason  for  suggest¬ 
ing  this  approach  is  quite  simple:  the  alternative  of  calculating  the  required 
validity  conditions  (those  in  Appendix  B  are  just  the  start)  would  be  unreason¬ 
ably  lengthy  and  complicated.  It  is  therefore  impractical  and  unnecessary  to 
present  a  list  of  validity  conditions.  The  justification  for  the  method  suggested 
lies  in  the  assumption  that  moderate  changes  in  parameters  from  the  nominal  values 
will  not  affect  the  segregation  of  large  and  small  terms  (except  for  Configura¬ 
tion  4) . 

C.  ADEQ.UACI  OF  ONE-  AND  TWO-ELASTIC-MODE  REPRESENTATIONS 

Regardless  of  the  validity  of  the  approximations,  there  i.»  rtill  a  basic 
question  as  to  the  number  of  modes  required  to  adequately  represent  the  system(s) 
under  study.  This  subject  was  investigated  in  Ref.  1  for  the  three  cases  treated 
here  and  the  results  of  this  investigation  are  summarized  below. 

Configurations  £  iu'.d  3  shown  to  be  accurately  represented  with  only  one 
or  two  flexible  modes,  the  frequency  response  curve  being  accurate  ^as  determined 
by  coD^arison  with  a  five-elastic-mode  case)  up  to  the  characteristic  frequency 
of  the  Ijist  flexible  mode  included.  However,  Configuration  4  was  shown  to  require 
considerably  more  flexible  modes  for  an  accurate  representation.  In  an  effort  to 
obtain  some  usable  data.  Configuration  4  was  investigated  at  several  conditions 
o.';  ro'^uced  dynamic  pressures  lc-''°r  than  those  studied  in  Ref.  1 .  These  lower 
preoc  reo  cei.a  <.o  minimize  the  dynamic  effects  of  the  higher-frequency  modes,  as 
indicated  by  the  small  changes  which  occurred  in  ’ch  ■  exact  factors  for  tlie 
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TABLE  II 

TRANSFER  FUNCTION  APPROXIMACE  FACTORS 
CONFIGURATION  2 


TABLE  III 

TRANSFER  FUNCTION  APPROXIMATR  FACTORS 
CONFIGURATION  3 


TABIE  IV 

TRANSFER  FUNCTION  APPROXIMATE  FACTORS 


TABLE  V 


TRAMSFEK  FUNCTION  APPROXIMATE  FACTORS 
CONFIGURATION  2 

4  MODES 
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TABLE  V  CONCLUDED 
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TABLE  VI 


TRAMSFER  FUNCTION  APPROXIMATE  FACTORS 
CONFIGURATION  3 

4  MOOES 
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’’’ABLE  VI  CONCLUDED 


transfer  FtJNGTiON  APPROXIMATE  FAC?rORS 
CONFIGURATION  4 
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three-  and  four-mode  transfer  functions.  The  range  of  dynamic  pressures  for 
which  the  one-  or  two-elastic-mode  representation  is  valid  (regardless  of  the 
validity  of  the  approximate  factors)  is  thus  strongly  limited  for  this  one  con¬ 
figuration.  Care  should  he  exercised  if  this  type  of  configuration  is  to  be 
represented  with  onlj'  a  few  of  its  normal  modes. 

D.  NUMERICAL  COtffARISONS  OF  EXACT  AND  APPROXIMATE  FACTORS 

The  excellent  agreement  between  the  approximate  and  the  exact  factors  of 
Configurations  2  and  3  indicates  that  the  approximation  formulas  for  these 
configxxrations  can  he  expected  to  remain  valid  for  extreme  ranges  in  dynamic 
pressure.  However,  the  approximation  formulas  for  Configuration  4  were  not 
valid  when  the  dynamic  pressure  was  extended  to  20  psl.  Thus,  those  approxima¬ 
tions  should  he  used  cautiously  when  conditions  of  high  dynamic  pressures  are 
investigated. 
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SECTION  IV 


SINGIE  SENSOR  CONTROL  LOOP  SYSTEMS 


A.  INTRODUCTION 

Control  of  the  longitudinal  axis  in  general  Implies  control  of  the  short 
period.  Accordingly,  the  closed- loop  bandwidth,  or  equivalently  the  open- loop 
crossover  frequency,  must  be  roughly  equal  to  or  greater  than  the  short-period 
frequency.  Vfhen  flexibility  effects  are  present  such  crossovers  can  sometimes 
lead  to  closed-loop  instabilities  because  of  structural  "coupling"  excited  by 
the  autopilot.  Such  incipient  instabilities  can  easily  be  investigated  by 
Bode  analysis,  and  can  in  general  be  avoided. 

Since  the  output  quantity  fed  back  to  the  controller  senses  all  conrponents 
of  motion,  rigid-body  as  well  as  elastic  (unless  filtered),  the  nature  of  the 
complete  open- loop  tmnsfer  function  can  often  be  drastically  changed  by  a  change 
in  sensor  location.  Thus,  whereas  for  a  given  sensor  location  it  may  be  impossible 
to  cross  over  near  the  short-period  frequency  without  appreciable  excitation  of  an 
elastic  mode,  a  slight  shift  in  sensor  location  may  permit  reasonable  closures. 

Both  the  general  formulation  of  the  output  quantity  as  a  function  ol’  sensor 
location  and  the  process  of  selecting  a  "proper"  location  are  discussed  below, 
with  specific  reference  to  the  use  of  vertical  gyro  feedback  loops. 

B.  SENSOR  OUTPUT 

If  the  four-mode  perturbation  equations  involve  the  vertical  displacement,  h, 
the  pitch  angle,  0,  and  the  first  two  elastic  modes,  and  gj,.,  then  the  air¬ 
frame  transfer  functions  are  h/6,  0/5,  gj/S,  and  gk/S.  The  rigid-body  degrees 
of  freedom  a^'e  h  and  0,  while  g^  and  gij  represent  the  first  two  elastic  degrees 
of  freedom. 

The  deflection  at  any  point  i  along  the  fuselage  reference  line  will  be 

yi  =  q)^,h  +  9120  +  (Pijlj  +  <Pi4l4  (75) 

T*!"  --ilope,  or  pitch  aa<Jl..  measured,  for  example,  by  a  vertical  gyro  at  any 


point  i  is  found  by  differentiation: 


%  = 


^'Pi^  ‘^‘Pfl^ 

dx  dx  dx  ^5  dx 

<Piih  +  'Pia®  -Pijb  ■" 


But  for  the  rigid-body  modes, 


“Pil  =  ° 
fia  -■=  ^ 


Thus, 


0i  =  e  + 


(76) 

(77) 

(78) 


Because  the  sensor  will  detect  the  total  physical  motion,  the  ransfer  function 
which  must  be  considered  is 


‘'i 


t  ■*■  *Pi5  ‘Pi4  B 


(79) 


■following  the  transfer  function  factored  forms  given  in  Table  I, 
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Therefore, 


fi 

5 


»e  ,  .  \ 

A  •^'Pi3-^  +  'Pi4-A- 


(81) 


For  convenience,  this  can  be  written  in  a  slightly  different  fom: 


0i  r  .  N^ 

'8  =  [  'Pi?  %  ®i4  N J  Z 
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(82) 


The  sensed  motion  is  thus  the  mean  centerline  motion  (0)  modified  by  the 
bracketed  factor  of  E(i  82.  The  mean  centerline  transfer  function  is,  itself, 
different  from  the  rigid-body  transfer  function  obtained  when  the  elastic  moles 
are  neglected.  The  difference  is  apparent  from  the  following  equations  (cf. 
Table  I). 

Rigid  Alrfrm." 


Elastic  Airframe  (Two  Elastic  Modes) 


3[s2  +  (25a3)5pS  +  a|p]p  +  (25a))^^s  +  [s^  +  (2ta>)2gS  +  a|e] 


(84) 


The  addition  of  two  elastic  modes  lias  added  two  pairs  of  second-order  roots  to 
the  numerator  and  denominator  of  the  mean  centerline  response,  0/5.  This  is  in 
addition  to  the  elastic  inputs  which  are  added  to  0/5  as  shown  in  Eq.  79* 

C.  CLOSED-LOOP  CONSIDERATIONS 

Assuming  lead  equalization  of  the  gyro  output  (T(,is  i  ')  and  neglecting 
actuator  and  sensor  lags,  the  open- loop  transfer  function  given  by  Eq  84  yields 
a  frequency  response  curve  of  the  form  shown  in  Fig.  4.  The  actual  curve  will 
depend  on  the  relative  positions  and  degree  of  damping  of  each  pair  of  complex 
roots  shown.  Closing  the  loop  in  this  case  is  quite  simple,  requiring  only  that 
the  zero  db  line  inlr-'^',..  .aapZ-itude  curve  at  a  frequency  in  the  neighbor¬ 
hood  of  <%p  and  yet  not  intersect  either  of  the  higher- frequency  "peaks."  (Servo 
lags,  unimportant  at  short  period,  will  reduce  the  phase  margin  at  higher  frcq"..-ii~ 
cies.)  Such  a  closure  would  be  impossible  if  the  0/5  frequency  response  curve 
were  of  the  form  shown  in  Fig.  5  where  a>]e  is  assumed  lower  than  (and  is  very 
lightly  damped) .  For  this  latter  case,  crossover  near  would  result  in 
i’>stub:l...  y  near  o, 
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Figure  4.  Bode  (jco)  Amplitude  for  Typical  e/6  with  Lead  Equalization 
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If  the  0/6  response  is  favorable,  as  shovm  in  Fig.  4,  then  the  complelie  0i/5 
response  can  also  be  made  favorable  by  minimizing  the  dynamics  of  the  bracketed 
tenn  in  Eq,  82  (via  the  sensor  location) .  If  the  0/8  response  is  similar  1:o  that 
shown  in  Fig.  5,  however,  then  the  sensor  should  be  located  so  the  dynamics  of 
the  bracketed  term  reshape  0i/S  to  obtain  a  more  desirable  form. 

D.  OPTIMUM  SEKwJR  LOCATION 

Methods  for  locating  the  sensor  to  achieve  desirable  results  for  the  two 
postulated  situations  will  now  be  described.  In  the  first  case  the  sum  of 
(PIjNij/Nq  and  <Pi4Nj^/Ng  will  be  held  approximately  constant  over  the  frequency 
range  of  interest.  Equation  82  shows  that  under  these  circumstances  0i/8  will 
be  equal  to  0/8  with  a  gain  change.  In  the  second  case  a  sensor  location  will 
be  found  which  makes  the  combined  dynamics  of  the  terms  in  the  brackets  of  Eq  82 
Just  cancel  the  elastic  modes  found  in  the  mean  centerline  transfer  function, 

0/8  (Eq  84) . 

The  example  chosen  to  demonstrate  these  methods  is  associated  with  the  high 
q  condition  (1197  psf)  for  Configuration  J.  The  transfer  functions  given  in 
Eq  82  can  be  obtained  from  Appendix  D  and  are  repeated  below.  (Note  that  the  0/8 
transfer  function  will  produce  a  frequency  response  similar  to  that  shown  in 
Fig.  4  and  therefore  is  satisfactory.) 

=  (s  +  1 .53)  [s^  +  2(0.25)(12.8)s  +  l65][s2  +  2(0.015)(26)s  +  6?^ 

=  s[82  +  2(0.22)  (15. 1  )s  +  228]  [s2  +  2(0.39)  (77.2) +  5987] 

To  -ir  -1  (®5) 

=  s|_s2  +  2(0.17)  (3.1 6)s  +  26.7j[82  +  2(0.23)(l5.l)5  +  17^ 

=  b[b2  +  ;,,u.:»2)i4.^y)s  ...  .1]  7(0.22) (l 2.9) s  +  l6^ 

X  [b2  +  2(0.081  )(29.7)s  +  884] 


A0 
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The  process  of  locating  the  sensor  to  make  the  bracketed  terms  of  Eq  82 
independent  of  frequency  is  simplified  by  considering  an  alternate  form. 


r 
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(86) 


From  Eq  86,  it  is  obvious  that  mthamatl  cal  operations  equivalent  to  closing  two 
loops  will  be  performed,  i.e.,  must  be  added  to  1 .0  and  the  result 

must  be  multiplied  by  <p^,.N^i. /Ng  and  added  to  1.0.  The  solid  line  in  Pig.  6 
shows  th>»  Bode  (Jod)  plot  for  the  amplitude  of  and  Fig.  7  shows 

the  Bode  (joj)  plot  for  the  amplitude  of  AgNjj^/Ajj^Ng.  When  the  gain  for  the 
first  closure  (tp]  A|j/(p|[^A»ij)  is  chosen,  then  the  curve  for  1  + 
will  resemble  that  shown  by  the  dashed  line  in  Fig.  6.  This  is  easily  seen  by 
noting  that  for  any  transfer  function,  G, 

1  +  G  A  G  for  G  »  1 

1  +  G  A  1  for  G  «  1 

The  closed- loop  curve  (corresponding  to  closing  the  first  loop),  is  therefore 
closely  approximated  by  when  that  quantity  is  much  greater  than 

unity  (zero  db),  and  by  the  zero  db  line  when  <Pi  is  much  less  than 

unity.  For  regions  where  G  A  i ,  the  closed-loop  can  most  conveniently  be  plotted 
using  conventional  Nichols  charts. 

Simultaneous  inspection  of  the  dashed  line  of  Fig.  6  and  the  plot  in  Fig.  ^ 
shows  that  the  two  curves  have  a  mirror  i-mage  resemblance.  Ibis  is  a  result 
of  a  Judicious  closure  of  the  first  loop  (i.e.,  properly  locating  the  zero  db 
line  in  Fig.  6) .  Because  these  two  curves  represent  quantities  which  are  to  be 
multiplied  (and  thus  their  logarithmic,  db,  plots  are  to  be  added),  the  product 
is  seen  to  be  relatively  independent  of  frequency,  r  e  3;ey  point  here  is  that 
the  gain  of  the  first  closure,  ,  was  chosen  to  appropriatelj'  locate 

the  zero  db  line  in  Fig.  6.  Tlie  corresponding  appropriate  sensor  location  can 
now  be  determined  from  Fig.  8,  which  gives  the  value  of  (p^jA|j/q)]^I^A||^  as  a  func¬ 
tion  of  fuselage  station.  Detailed  numerical  considerations  show  that  the 
sensor  should  be  placed  station  656  in  order  to  make  the  value  of 


relatively  independent  of  frequency.  Closing  the  second  loop  is  now  trivial 
ba;»usc  j.t  ..>orr"cpo:id':  adlv'.  a  constant  to  1 .0. 
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It  is  noted  that  the  mirror- Image  effect  requires  that  a  zero  of  N|.^  /Ne  and 
a  pole  of  occur  at  the  same  frequency  with  the  same  damping.  However, 

this  is  always  the  case,  because  the  zeros  of  the  former  are  the  poles  of  the 
latter.  This  method  will  therefore  always  theoretically  reduce  Si/Bg  to  approx¬ 
imately  v3/6g  (when  only  two  elastic  modes  are  considered). 

For  the  second  situation  postulated,  i.e.,  a  @/6  frequency  response  as  shown 
in  Fig.  5<  the  dynamics  of  the  elastic  modes  can  be  used  to  better  shape  the 
sensed  pitch  response.  (Note  that  no  such  modification  is  required  for  the 
example  picked.)  Again,  there  are  two  closures  involved  (see  Eq  86),  with  the 
location  of  the  zeros  for  the  second  closure  depending  on  the  gain  assoclatea 
with  the  first  closure.  In  turn,  the  gain  is  strictly  dependent  on  the  sensor 
location  (mode  shapes). 

If  the  zeros  resulting  from  the  two  closures  are  to  be  placed  in  close 
proximity  to  the  elastic  roots  in  A  (Eq  85)  with  a  resulting  cancellation,  the 
following  considerations  apply:  The  roots  of  the  first  closure  will  be  the 
zeros  for  the  second  closure;  one  pair  of  these  will  be  lightly  damped  and  close 
to  13  rad/sec  for  any  value  of  gain,  as  may  be  seen  by  inspecting  Fig.  9*  Since 
the  roots  of  Hg  also  Include  a  pair  near  13  rad/sec,  the  final  zeros  (which-are 
the  roots  of  the  second  closure)  V'ill  include  a  lightly  damped  pal  r  at  approxi¬ 
mately  13  rad/sec  (see  Fig.  10).  This  is  true  because  a  pole  and  zerc  in  close 
proximity  will  always  yield  a  i»ot  in  that  neighborhood  for  all  values  of  gain 
(provided  the  remaining  poles  and  zeros  are  relatively  far  removed,  as  they  are 
in  this  case).  Thus,  for  the  example  chosen  the  elastic  poles  of  A  at  approxi¬ 
mately  13  rad/sec  will  be  cancelled  for  any  sensor  posi.-cion  selected.  The 
placement  of  the  sensor  can  thus  be  made  with  the  intention  of  producing  a  pair 
of  lightly  daaped  zeros  at  approximately  the  location  of  the  second  elastic 
poles  of  A,  29.7  rad/sec.  The  second  closure  has  a  pair  of  lightly  damped  poles 
at  approximately  26  rad/sec.  The  locus  of  the  roots  emanating  from  these  poles 
must  depart  in  thr  In  PI5  10  if  the  locus  is  to  include  the 

desired  location  (29.7  rad/sec)  for  the  zeros.  Th;  poles  of  che  second  closure 
(being  the  Ng  numerator)  are  not  a  function  of  sensor  location;  and  of  the  four 
zeros  resulting  from  the  first  closure,  two  are  essentially  independent  of  sensor 
location  (the  two  lightly  damped  roots  at  approximately  13  rad/sec).  Thus  the 
problem  is  reduced  to  closing  the  first  loop  so  that  among  the  resulting  loo-.-s 
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(which  will  be  zeros  for  the  second  closure)  there  will  be  a  complex  pair  of 
zeros  which  yield  the  proper  second-closure  departure  from  the  first-closure 
poles  at  26  rad/sec.  Figure  10  shows  the  required  location  of  these  zeros 
(cross-hatched  area).  There  are  two  possible  ways  of  producing  roots  in  the 
croF,8-hatched  area:  the  gain  of  the  first  closure,  must  be  either 

positive  or  an  extremely  large  negative  number.  The  root  locus  for  the  first 
closure  is  sh—Ti  in  Fig.  9.  If  the  gain  is  positive,  it  must  be  large  enough 
to  yield  roots  greater  than  29.7  rad/sec.  With  this  information,  a  number  of 
positive  gains  are  tried  until  the  zeros  resulting  from  the  first  closure  fall 
in  the  desired  position  for  the  second  closure.  For  the  exanple  chosen,  the 
roots  of  the  first  closure  need  to  be  driven  all  the  way  to  the  zeros.  There¬ 
fore,  the  sensor  should  be  located  at  station  Gk'J  where  the  gain  of  the  first 
closure  takes  on  its  largest  value. 

It  is  unfortunate  that,  for  the  example  chosen,  both  methods  result  in 
placing  the  sensor  at  a  fuselage  station  where  the  gain  of  the  first  closure, 

is  changing  quite  rapidly.  Any  departu.e  of  the  mode  shapes  from 
those  expected  may  not  result  in  the  required  gain  and  hence  not  achieve  the 
desired  flat  response  or  the  required  cancellation  of  roots.  The  effects  of 
such  perturbations  on  the  system  can  easily  be  estimated  by  applying  the  Bode 
techniques  of  Fig.  6  and  7*  In- any  event,  regardless  of  the  validity  of  the 
exauples  chosen  to  illustrate  the  two  approaches,  the  basic  considerations 
involved  are  generally  applicable  to  the  closed- loop  analysis  and  synthesis  of 
flight  control  systems  for  flexible  airframes.  It  is  especially  pertinent  to 
note  that  a  fairly  complete  set  of  transfer  func‘‘icns  (.ln''luding  those  of  the 
coxpled  elastic  modes  themselves)  is  required  for  such  analysis  and  synthesis 
activities. 
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SECTION  V 


RECOMMENDATIONS  FOR  FUTURE  INVESTIGATIONS 


The  Gtudy  leading  to  the  results  presented  in  this  report  has  disclosed,  or 
has  investigated  only  partially',  several  potential  areas  of  research,  including 
the  fcllotfing: 

1 .  the  prohlcm  of  finding  the  quasi-static  aeroelastic  correction 
for  equations  of  motion  which  include  a  few  elastic  modes 

2.  the  basic  mechanics  of  mode  interaction,  and  an  understanding 

of  what  parameters  can  best  be  expected  to  provide  an  indication 
of  the  degree  of  coupling  present 

5.  the  possibility  of  representing  the  motion  of  an  elastic  air¬ 
frame  with  a  simplified  set  of  equations  of  motion 

if  the  transfer  function  approximate  factors  are  a  lanction 
of  mode  shape,  what  approximations  can  be  made  to  adequately 
approximate  the  required  modes? 

5*  when  and  how  can  the  approximations  developed  in  this  study 
best  be  utilized? 

As  indicated  by  the  manipulations  outlined  and  demonstrated  in  Appendix  A, 
any  physical  feeling  for  the  meaning  of  [Xcq]  is  con^letely  lost  in  the  maze 
of  relationships  involved.  This  places  the  method  in  the  category  of  being 
completely  unsuited  for  use  in  the  practical  calculation  of  approximate  transfer 
function  factors.  Nevertheless,  the  basic  feeling  e^sts  that  in  ,.ie  way  or 
another,  [x®]  must  cci'rcspor.d  to  "modified  elastic  properties."  That  is, 
the  number  of  elastic  modes  already  included  must  give  rise  to  a  correction 
of  the  basic  static  deflection  properties;  and  the  corrected  properties  (the 
"residual  stiffness")  must  in  some  way  be  connected  with  [Xqo]  •  These  con¬ 
nections  have  to  be  for.r.ed,  and  a  physically  sutisfyiog  picture  must  be 
drawn  before  the  process  of  obtaining  simplified  approximations  to  xhe  aero- 
olastic  corrections  can  proceed. 

The  desirability  of  obtaining  such  approximate  corrections  cannot  be  over¬ 
emphasized.  This  stems  basically  from  the  fact  that  if  procedures  akin  to  those 
described  in  Section  II  are  required  to  establish  proper  aeroelastic  corrections, 
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then  much  of  the  iini)act  of  approjcimate  transfer  functions  is  lost,  hecause 

1 .  the  time  and  the  machine  methods  required  to  compute  [X  ool 
might  just  as  well  he  used  to  compute  exact  transfer  functions 

2.  the  basic  possible  physical  insight  which  approximate  transfer 
functions  can  potentially  yield  will  not  have  been  realized. 

A  second  major  consideration  is  the  phenomenon  of  coupling.  Better  under¬ 
standing  of  coupling  is  required;  as  evidenced  by  the  results  presented  here 
on  Configuration  5.  In  that  instance,  the  equations  proved  to  be  weakly  coupled, 
even  though  it  was  expected  that  strong  coupling;  would  be  present.  The  reason 
for  this  is  now  know: :  coupling  of  two  modes  is  not  necessarily  indicated  by 
proximity  of  their  frequencies.  This  is  discussed  in  Ref.  8. 

The  equations  of  motion  for  flexible  vehicles,  including  two  rigid-body  and 
an  arbitrary  number  of  flexible  degrees  of  freedom,  arc  given  in  their  most 
compact  form  by  Eq  jk. 

The  approximate  transfer  functions  derived  in  this  study  are  a  direct  in¬ 
dication,  for  the  cases  studied,  of  the  relative  importance  of  the  various  terms 
in  Eq  74.  Unfortunately,  as  shown  in  Section  III,  all  terms  appear  in  one  or 
another  of  the  various  factors  involved  in  the  complete  set  of  transfer  functions 
Accordingly,  the  specification  of  the  validity  conditions  for  which  the  approxima 
tions  apply  becomes  exeedingly  complicated.  Because  all  the  parameters  remain 
important  (depending  on  the  particular  root  involved),  the  most  efficient  way, 
currently,  of  determining  the  applicability  of  the  approximations  is  to  conq)ute 
an  "exact"  check  case.  If  the  approximations  are  valid  a.-  snown  by  this  com¬ 
parison,  they  can  be  applied  to  gain  the  desired  insight  into  sources  of  dif¬ 
ficulty,  effects  of  changes,  etc. 

An  alternative  approach  to  approximate  factorization  is  to  write  sets  of 
"simplified  equaulo’ia  "o+ion, "  each  set  applicable  to  restricted-frequency 
regions  (e.g..  Kef.  9)*  The  sets  of  simplified  equations  and  ^he  sets  of 
approximate  factors  are  complementary  ways  of  specifying  the  important  conti'?'.,- 
uting  terms;  both  approaches  will  theoretically  yield  similar  results  for  the 
approximate  transfer  functions.  Accordingly,  the  simplified  equations  can  be 
used  to  specify,  hopefi^y,  more  tractable  validity  conditions,  and  fuither  t-, 
approximate  transfer  iXinctions  for  other  than  control  inputs  (e.g., 

.'oi  gusts,  as  in  Ref.  9).  For  these  reasons,  the  simplified  equations  are 
highly  desirable  adjuncts  to  the  appi'oximate  +-ansfer  function  factors. 
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The  approximate  factors  for  the  transfer  functions  of  a  flexible  airframe 
are  given  in  terms  of  such  quantities  as  airp)lane  stability  derivatives  and 
mode  deflections  at  various  points  on  the  airframe.  A  great  many  reports  and 
papers  deal  with  the  subject  of  approximating  such  parameters,  and  many  of 
those  concerned  with  stability  derivatives  are  directly  applicable.  The  same 
is  not  true  of  all  methods  of  estimating  mode  deflections,  most  of  which  are 
used  only  to  establish  an  initial  estimate  for  use  in  iterating  to  the  exact 
value.  Because  the  iteration  procedures  rapidly  converge,  the  initial  estimates 
need  not  be,  and  are  not,  very  accurate.  This  is  especially  true  of  Galerkin's 
iteration  method,  and  the  method  of  Stodola  and  Vianello,  where  the  iteration  is 
continued  until  repeated  iterations  provide  the  same  answer  (Ref.  5) .  Other 
methods,  such  as  Rayleigh-Rltz,  modified  Rayleigh-Rltz,  collocation,  and  colloca¬ 
tion  using  station  functions  (Ref.  3),  rely  somewhat  more  heavily  on  the  original 
estimate  if  any  accuracy  is  to  be  obtained.  Thus,  many  techniques  of  varying 
accuracy  (directed  at  these  latter  methods)  have  been  formulate  ',  to  provide  a 
fairly  reasonable  estimate  of  the  mode  shape.  Most  of  these,  however,  are 
usually  content  with  merely  satisfying  boundary  conditions.  The  importance  of 
selecting  or  developing  such  approximations  stems  from  the  fact  that  fairly  large 
errors  in  mode  shape  may  be  tolerable  for  the  purpose  of  computing  approximate 
transfer  function  factors .  The  effect  of  mode  shape  error  on  the  accuracy  of 
the  approximate  transfer  functions  is  easily  determined  in  a  given  case  by 
findir^  the  changes  in  the  factors  caused  by  variations  in  the  values  of  the 
mode  deflections. 

It  is  desirable  to  epply  the  approximate  transfer  function  fOiinulas  (Ref.  5) 
to  some  actual  aircraft  or  missiles  to  demonstrate  their  application  and  utility 
on  a  tangible  basis.  Probably  the  most  significant  and  useful  results  can  be 
obtained  for  vehicles  that  are  currently  in  the  preliminary  design  stage.  The 
information  necessary  for  !.!.<;  c'^luation  of  the  approximate  factors  is  generally 
most  available  at  that  time,  and  the  resulting  analysis  woidd  be  useful  to  the 
vehicle  manufacturer,  inasmuch  as  some  insight  would  be  provided  on  the  dominant 
factors  affecting  the  aircraft  modes. 

Many  of  these  problems  nay  be  resolved  in  the  completion  of  Contract  No. 

AF  33(657) -837^  which  has  recently  been  awarded  specifically  for  study  in  these 
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APPEKDIX  A 


AEROELASTIC  CORRECTIONS 


A  method  has  been  derived  in  Ref.  1  wherein  IXqq]  can  be  calculated.  It  was 
not  considered  necessary  to  repeat  that  derivation  in  this  report,  but  a  summary 
of  the  method  and  an  example  of  its  use  follow. 

A.  SUMMARY  OF  THE  METHOD 


To  obtain  the  basic  data  required  in  these  equations,  the  system  is  restrained 
at  two  points,  and  is  then  subjected  to  a  unit  acceleration  field,  first  in  trans¬ 
lation  and  then  in  rotation.  Qhe  resulting  physical  deflections  are  partitioned 
in  matrix  form  as 


•Jg  trans  I 


%  trans 


rot 
^  rot 


(A-1 ) 


(A-2) 


where  the  subscripts  a,  g,  and  m  denote  acceleration  deflections,  grounded  coordi¬ 
nates,  and  movable  coordinates,  respectively.  These  ‘'.eficctionr.  are  used  in  the 
calculation  of  rwi>  where 


In  Eq  A-5,  the  columns  of  correspond  to  the  finite -frequency  mode  shapes 
(those  of  interest  as  mentioned  earlier);  represents  the  diagonal  [k  +  Ms^ 
matrix  for  the  corresponding  finite  frequency  stiffness  and  mass  matrices;  and 
the  elements  of  [Xp+Q^]  are  found  by  calculating  partitions  of  [Xf4^]  as  follows: 


^(f4<o)gg  :  ^(f+co)gra 

D<fU  . ; .  (“-‘‘i 

4 

(f+C!o)ng  ^(f+cc)Km, 

[Hf^)4  ■  transit  rot}]  (^-5) 

|^*(f+QD)mg]  “  "^1*30  trans}|%  rot}]  [^qj] 

K]  [Mo]  ■’  [^of[^  [{<3a  trans}  rot}]  [“o]  [%]''  (A-6) 

j^^(f+ao)gm]  “  |^^(f+oo  )nig]  (A-7) 

[^(f+aD)inm]  ”  |^^]  *  [*°®](5og]  [^(r+<»  )gg]  [*og]  [*0“] 

^  [^(f+CD)gm]  ■•■  [^(f+oo  )ing]  ['*’og]  ['*’om] 


where  C2b]  is  the  influence  coefficient  matrix  of  the  cysrem  when  restrained  at 
the  two  points  qg)  and  where  jjj*oin]  ail'll  [^^s]  found  by  partitioning  the  zero- 
frequency  modal  matrix  into  the  elements  corresponding  to  the  restrained  and 
unrestrained  coordinates, 

a] 

Also,  is  the  zero-frequency  modal  mass  matrix 

H  is  the  mass  matrix  of  the  physical  system  (defined  previously) 
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B.  EXEMPLARY  DEVELOPMENT 

The  three  masses  shown  in  Pig.  A-1  are  rigidly  attached  to  a  weightless  beam 
of  length  21: 


Figure  A-1 .  Mechanical  Model 


The  rotational  inertia  of  the  first  and  third  masses  in  Fig.  A-1  is  negligible, 
and  the  rotational  inertia  of  the  second  mass  is  Ig*  All  deflections  are  measured 
inertially,  and  are  positive  as  indicated.  The  angle  0  is  the  inertial  rotation 
of  mass  two,  and  is  also  positive  as  shown.  Mass  two  is  in  the  center  of  the  beam 
which  has  a  constant  El.  The  masses  are  assumed  to  be  equal,  m. 

With  this  information,  it  will  be  possible  to  write  the  system  equations  of 
motion  utilizing  Lagrange's  equation: 


d  /  5l\  dL  ^  5D 

dt  WrI  Wr 


Qj. 


(A-10) 


where 


L  =  T  -  U 

T  =  the  kin:,  vie  energy  of  ll:  •system 
U  =  the  potential  energy  of  the  system 
D  =  the  damping  energy  of  the  system 
Qj.  =  the  generalized  force  on  the  r'*'^  degree  of  freedom 
q^  =  the  r^h  degree  of  freedom 


Follow ‘.;ig  rufme  drfi:;it!o:«tj  (an  session  for  the  potential  energy  Tor  the  beams 
is  derived  in  Ref.  J) , 
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L  =  I  +  I  BiJg  +  I  “yf  I  ^2^ 

- 1  ^  (y,  -  yg  -  ie)2 

- 1^  (yj  -  y2  -  10)^  (A-11) 


It  is  not  necessary  to  consider  any  external  forces  for  the  purposes  of  this 
exan5)le;  thus,  Qj.  will  be  zero.  Substituting  in  Lagrange's  equation  yields 


®y2  ■  5  Q  (yi  -  y2  -  10)  -  3  ^  (yj  -  y2  +  10)  =  0 

igO  -  3  P  (y,  -  y2  -  10)  +  5  ^  (yj  -  y2  -  10)  =  0 

FT 

wj  +  3  ^  (y3  -  y2  +  10)  =  0 

ffij'i  +  3  ^  (y^  -  y2  -  10)  ==  0 


(A-12) 


Equations  A-12  can  be  put  into  matrix  form  as  indicated  in  Section  II: 

[ms2  +  k]{qf  =  Iq}  (A-13) 


where 


0 

0 

0 

m 
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A 

\J 

'  l3 

'l3 

r  El 

*  El 

X  El 

[k]  = 

0 

b-5- 

'32 

-'32 

-  El 

,  El 

»  El 

■•'13 

53 

3  “x 
l3 

0 

,  El 

X  El 

[■'13 

-3  -0 

1^ 

0 

'13 

The  modal  matrix  can  now  be  found: 

1 

0 

-2 

0 

0 

1 

0 

-2?i 

1 

I2 

1 

-1 

1 

-1 

.  1 

1 

1 

1 

(A-14) 


If  the  degrees  of  freedom  of  Interest  were  limited  to,  for  example,  the  two  rigid- 
body  modes  and  just  one  elastic  mode,  then,  from  Section  II, 


1  0 

°  I 

1  -1 


1  1 


r.,  1  _ 

I  'T  • 

_  t 


-2 

0 


(A- 15) 
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The  [y]|  matrix  can  also  te  calculated  ao  described  in  Section  II,  and 


Following  Section  II, 

*  308^  0 

H  =  0  (2in  +  ^)s2 

[Yf]  =  [^6ms2+5lf||j 


(A-17) 


From  this  information,  j^^ooj  calculated  direcoi^  fiom  Eq  47: 

[x„]  .  WH'Jv]*  ("7) 


(A-18) 


'^1 


As  described  in  A-1  through  A-<),  jX^^J  con  also  be  calculated  without  the 
knowledge  of 

M' 

With  the  system  restrained  in  translation,  and  rotation  at  the  center  of 
gravity  (mass  two) : 


6a 


Equation  A- 6  yields 
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From  Eq  A-7, 


27EI 

2i^l^  +  2nUJ'2 

4i5 

2m2l5  +  2mllp 

.L  ^ 

3Si(a.  ♦ 

201^15  +  2ralIo 

27EI 

|^(f+oo  )dbi^  “ 

4i5 

3El(2m  + 

2m®l5  +  Emllg 

“  '^TEI 

2 

5El(2m  +  ^ 

27EI 

3El(2.  + 

Now  [X(f+oo) 

can  be  constmc 

ted  from  Eq  A-21 ,  A-23, 

A-24,  and  A-26j 

[*(f+co]  = 

21^ 

l3 

l5 

27EI 

0 

■  27EI 

■  27EI 

0 

Em^l 

mig 

mIg 

5Ei(a.  *  ^.) 

I  ^ 

5El(2m  + 

lo  ^ 

3El(2m  +  J) 

l3 

mig 

5j^3  2m^l3  +  2mlTg 

1^2_3  2m^l5  +  2mll2 

■  27EI 

j  2 

3EI^2m  + 

27EI  2 

3EI^2m  + 

27EI  '  I  .2 

3Erf2m  + 

l3 

ml2 

1j]^3  2m^l3  +  2mll2 

5^3  Pm^l^  +  2mll2 

■  27EI 

2 

■  2TEI  '  I  2 

27EI  .  I.  2 

5El(2m  +  ^) 

3EI(2d.  +  ~j 

3El(2m  +  -|) 

\  1^' 

Before  solving  for 


from  Eq  A-5^  it  will  be  necessary  to  find 


from 


(A-28) 
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l5 

55ei 


(A-29) 


Subtracting  Eq  A-29  from  Eq  A-27, 


k  0  -2  -2 

0  0  0  0 

-2  0  1  1 

.-2  0  1  1  . 


0 

0 

0 

0 

2m2l 

mig 

ml2 

KJ 

2 

5 

/  lp\ 

3EI^2m  +  ^ 

1 

3El(2m  +  -5) 

3EI^2m  + 

fx  1  - 

A 

ml2 

l3  2m2l5  +  2mll2 

_  j2. 

2ni2l5  +  2mll2 

[XcoJ  - 

U 

/  Ip' 

"2 

^  ’  /  l2\^ 

"  6EI 

/  l2\^ 

3Eipm  +  ^ 

) 

3El(2m  + 

3El(2m  - 

A 

ml2 

l5  2m2l3  +  amllp 

i3 

2m2l3  +  2mll2 

U 

/  ^2' 

2 

6EI 

^  2 

3El|2m  +  ^ 

I 

3El(2m  + 

*  55) 

(A-30) 

Algebraic  manipulation  yields 


'0 

0 

_  0 

:  0 


0 

0 

0 

2m2l 

raig 

mig 

3EI^2m  + 

--f 

l2/ 

jEi  (a.  + 

5Ei(a.  * 

inl2 

I2 

li 

if 

i"/ 

j  2 

Aiw-i/gm  +  -7) 

/  ^2\^ 
6IEI |2m  +  p j 

raig 

li 

li 

3EI^2m  + 

-f 

l2/ 

I  ^ 

6lEl|2m  +  -|  j 

l2  ^ 

6lEI^2m  +  -|) 

(A- 31) 
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Equation  A-JI  reduces  to 


The  result  obtained  in  Eq  A- 52  proves  to  be  identical  to  the  result  obtain  in 
Eq  A-l8,  although  no  use  has  been  made  of  the  fourth-mode  shape. 


APPENDIX  B 


ANALYTICAL  METHODS  OF  APPROXIMATE  FACTORIZATION 


A.  CHARACTERISTIC  POLYNOMIAIB  OF  LIGHTLY  COUPLED  SYSTEMS 


Approximate  factors  of  a  characteristic  polynomial  are  found  directly  from 
the  matrix  of  coefficients  in  the  equations  of  motion,  rather  than  hy  expanding 
the  determinant  of  coefficients,  and  then  factoring  the  resulting  polynomial. 

The  technique  employed  here  involves  determining  those  corrections  that  must  he 
applied  to  a  crude  first -approximation  to  the  factors.  This  method  is  particu¬ 
larly  suitable  to  those  cases  where  the  diagonal  elements  of  the  determinant  of 
coefficients  are  the  major  contributors  to  the  characteristic  polynomial.  The 
case  to  be  considered  here  is  that  of  a  3-degree-of- freedom  system  having  a  deter¬ 
minant  of  coefficients  of  the  following  form: 


s  +  a„ 

®12 

Hf  ®13 

CO 

III 

®21 

s  +  agg 

^23^  +  ®’23 

®31 

^32 

q2  'bjjs  + 

(B-1) 


where  A(s)  is  the  characteristic  polynomial,  and  a^^j  and  bj^j  are 
real  constants 

This  could  represent  an  airframe  with  two  rigid-body  degrees  of  freedom  and  one 
elastic  structural  mode  (as  per  Eq  74).  A(s)  can  be  e:q)aii.  cd  about  the  third 
column,  giving 


where  p(s) 

Rl(s) 

and  ^2^®^ 


A(s)  =  P(s)  +  R^'s)  +  RgCs) 


+  by^.n  +  a:^j) 


s  + 


ai2 


■(^23®  +  ®23^ 

(b^js  +  aij) 


s  + 

^51 

^31 


s  +  agj 

®32| 


s  + 


^22 


a.'zo 


(B-2) 
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P(s)  can  also  1)6  written 


P(s)  =  (s^  +  Ps  +  y^)(8^  +  as  +  x^) 

=  8**'  +  (a  +  p)s^  +  (x^  +  .Y^  +  03)8^  +  (oy^  +  Px2)s  +  (B-3) 

=  s^  +  Bs''  +  Cs^  +  Ds  +  E 


where  by  direct  con^arison 

B  =  a  +  P 

C  =  x^  +  +  ap 

D  =  ay^  +  Px^ 

E  =  x^y^ 

and  by  comparison  with  Eq  B-2, 

x^  =  a^^a22  -  ^12^21 

,  (B-5) 

y-  =  ajj 

P  =>  bjj 

The  complete  A(s)  is  also  of  the  form  of  Eq  B-5,  but  with  sH.^htly  modified  fac¬ 
tors  and  polynomial  coefficients  due  to  the  added  R';(s)  and  R2(s/«  necause  the 
modified  polynomial  coefficients  are  directly  avails  ^le  from  Sq  B-2,  it  is  perti¬ 
nent  to  relate  increments  in  the  coefficients  to  increments  in  the  factors.  Then 
the  approximate  factors  of  A(s)  will  be  the  factors  of  P(s)  as  modified  by  these 
increments.  Proceedintj  al.yi.g  <.‘''*se  lines  by  taking  differentials  in  Eq  B-4, 


dB  =  da  +  dp 
dC  =  dx?  +  dy^  +  adp  +  pda 
dD  =  ody^  +  y^da  +  pox'”^  +  x^p 
dE  =  t  y^dx^ 


(B-6) 
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Because  R^(8)  +  R2(s)  jlS  of  second  decvee,  dB  =  0,  vrtiereBy 


da  =  -dp 


(B-7) 


Substituting  Eq  B-7  into  Eq  B-6, 


dC  =  dx^  +  dy^  +  (a  -  p)dp 
dl»  =  ady^  +  pdx^  +  (x?  -  y^)dp 
dE  =  x^dy^  +  y^dx^ 


(B-8) 

(B-9) 

(B-10) 


Eliminating  dp  by  combining  Eq  B-8  and  B-9  gives 

dC  -  dx^  -  dy^  dP  -  ady^  -  Pdx^ 
a  -  P  _  y2 


(B-n) 


Solving  Eq  B-10  for  dx  , 


d:^  =  ^ 


Substituting  Eq  B-1 2  into  Eq  B-1 1 , 


(B-12) 


dC 


-  <1,2  an  -  <*2  -  u 


a  -  P 


x"  -  •/ 


(B-1 3) 


Then,  solving  Eq  B-1 3  for  dy^, 


^,,2  =  (y^  -  :^)(y^dC  -  dE)  +  (g  -  P)(y^dD  -  pdE) 

+  (a 


(B-1 4) 


By  solving  Eq  B-10  for  dy^,  and  by  substituting  the  result  into  Eq  B-11,  a  E.’.cilar 
e:q)ression  for  d:^  can  be  found: 

^  _  (y^  -  x^)(dE  -  x^dC)  +  (g  -  p)(adE  -  x^dP)  {B-’3) 

(y^  -  x^)^  +  (g  -  P)(gy^  -  px?) 


71 


Expressing  Eq,  B-l4  and  B-15  in  terms  of  finite  differentials  (rather  than  infin¬ 
itesimal.),  and  collecting  teims. 


.2  _ 


[y2  -  +  a(a  -  -  [x?(y^  -  x^)]  ^  -  [x^(a  -  P)]ad 

(y^  -  +  (a  -  P)(ay^  -  Px?) 


(B-16) 


= 


•[y^  -  x^  +  p(a  -  P)]ae  +  [y^(y^  -  x^)]ac  +  [y^(a  -  p)]ad 


(y2  -  x^)  +  (a  -  P)(ay2  -  px2) 


(B-17) 


Considering  that  Eq  B-1  is  representative  of  an  airframe  with  one  elastic 
mode  include’d  in  the  equations  of  motion,  A(s),  in  terms  of  the  factored  forms 
of  Table  I,  would  be  given  by 


A(s)  = 


where 


(2yC)gpG 

'•  “tp]i 

s2  + 

sp 

=  x2 

+  Ax? 

(2?cs)sp 

^  w  H 

K  Aa 

2 

"^e 

-  y2 

+  Aj? 

(2!:a))ic 

=  P  ^ 

H  Ap 

(B-18) 


Therefore,  using  Eq-B-5,  B-l6,  and  .B-IT,  the  approximations  to  the  characteristic 
frequencies  are  given  in  terms  of  the  matrix  elements  by 

+  ®12^21  ®22Haii+a22-'b35)jAE 


■[(®11®22  ■  "  ®11®22  ®12'^2l)]^ 


I 

2 

%  ®'11®22  -  ®12®21 

i  -[(ann22  "  ®12^2l)(^11 +a22-b55)]AD  / 

^  (ajj  -  aiia22  +  \ 

i  ( 

(a,  .,  +  a22  -  i’55)[(®11  +  ^22^^53  '  ^55^®11®22  ‘  ®'12®21  ]  | 

(B-1 9) 
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0)^6  =  ajj-' 


I” [^^5  "  ^11^2  ®12®ai  ■*■  ‘’33(^11  °-22  ~  ^33^'^) 

'’’[^55^^53  ■  ^11^22  ■*■  a-|2a2l3-‘^  | 

+1^53(0.11  +  022  -  ^533^  ) 

,__  _ 

(055  -  0^^822  +  a,2a2i) 

I 

+(0^,  +  022  -  +  a22)a55  -  1555(011822  "  ai2a2i3 


r  (B-20) 


AC  =  -051I515  -  a^2^23 

AD  =  -0^j0j1  -  a25052  +  a5i(0i2l525  -  a22l5i5)  +  a52(a2il5i5  -  aiib25)  (B-2l) 
AE  =  aij(02i0j2  "  a22^31^  "  ®23(®’n^32  "  ®'12®’3l  ) 


Eq.u0tions  B-19  0nd  B-20  represent  first-order  corrected  values  for  the  squares  of 
the  short-period  end  first  elsstic-mode  frequencies  when  Bi(s)  end  R2(s)  ere  edded 
to  P(s)  to  give  A(s)  (Eq  B-2) .  However,  because  Eq  B-19  end  B-20  ere  very  un¬ 
wieldy,  it  is  desirable  to  simplify  the  two  corrections.  Subject  to  a  reesonsble 
set  of  velidity  conditions,  some  reletively  simple  relations  can  be  found.  Con¬ 
sider  the  following; 

Dividing  Eq  B-l4  by  Eq  B-15  gives 


-(y^  -  x^)(dE  -  y^dC)  -  (g  -  p)(pd^  -  :.^dD) 
(y^  -  :^)(dE  -  :^dC)  +  (a  -  S)(otdE  -  :^dD) 


(B-22) 


Dividing  numerator  and  denominator  by  (y^  -  x?)(dE  -  :^dC). 


This  can  be  greatly  simplified  by  making  the  following  assumptions  (which  ha:., 
nti  o5''’e’^'ed  to  be  time  ir.  many  instances): 

y2dC  «  dE 
x^dC  «  dE 
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Ktjuation  B-23  then  becomes 
-dx2 


(B-24) 


The  following  assumptions  have  also  proven  to  be  quite  reasonable  and  will  further 
reduce  the  complexity  "f  the  approximation: 


Further,  assurae 


whereby 


dx^  =  -dyS 


Substituting  Eq  B-25  into  Eq  B-10  gives 


(B-25) 


(B.  26) 


“sp  “  ®11®22  ■  ®12®’21 


ajj  -  an  822  + 


2  i  ,6E 

“^e  =  ®55  ■  a„  -  a^ao. 


^■33  "  1  ®22  -'2*^21 


(B-27) 


(B-28) 


Equations  B-27  and  B-28  are  the  desired  sinplifications  of  Eq  B-19  and  B-20. 

The  remaining  task  is  to  find  corrections  to  the  damping  terms  in  Eq  B-3. 
Because  the  frequency  corrections  are  now  known,  either  Eq  B-8  or  Eq  B-5  can  be 
solved  directly  for  the  tv  ons  which  (Sq  B-7)  are  simply  of 

opposite  sign.  Because  of  the  relative  magnitudes  of  the  quantities  .'nvclved, 
it  is  presumed  that  Eq  B-9  will  give  a  more  accurate  result  than  will  Eq  B-8; 
accordingly , 


dp  =  -da  = 


+  Pdx^ 


y^  -  xi^ 


(B-29) 


Traiisft  .'!."‘.ng  back  into  the  variables  of  interest,  and  replacing  differentials  with 
finite  differences. 


7k 


(^5^)gp  =  a-j  ^ 


AD  -  +  a22)Ay^  -  bjjAx^ 


22  *^35  -^1 1^22  '‘12*21 


(B-30) 


(25“)ie  =  1>53  - 


AD  -  (a^i  +  a22)Ay^  - 
*33  ■  *11*22  ■*■  *12*21 


(B-31) 


where  AD  io  given  in  Eq.  B-21  and  Ax?  and  Ay^  are  the  correction  terms  in  Eq  B-1 9 
and  B-20  or  in  B-27  and  15-28. 

It  is  now  desirable  to  have  a  simpler  expression  for  the  damping  con'ection. 
Using  Eq  B-25  in  Eq  B-29  gives 

dp  =  .te  4  (B.52) 


Now  Eq  B-26  can  be  substituted  into  Eq  B-32,  giving 


(B-35) 


Therefore,  Eq  B-30  and  B-31  cm  be  written 

AD  +  (; 


a^^  +  a, 


(25<s) 


sp 


=  a,  1  + 


11  *22 


V*53  -  ^1 


^22  '  °33  \ 

1®’22  *12^21' 


*55  "  *11*22  '’!2*21 


and 


AD  W-^1^1!221255 _ 

\^55  ~  ^11^22  *12^21  / 

®55  ■  ®11®22  ®12®'21 


(2Mie  =  ’’55 


(B-54) 


(B-55) 


The  following  is  a  summary  of  validity  conditions  which  allow  use  of  the 
simple  approximations  given  by  Eq  B-27,  B-28,  B-34,  and  B-55.  If  these  validity 
conditions  are  not  satisfied,  then  Eq  B-1 9,  B-20,  B-30,  and  B-31  must  be  used; 
in  this  case  only  the  first  validity  condition  listed  below  is  necessary. 
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Validity  Conditions 


1.  The  correction  terms  ai’e  all  small  (on  a  percent  liasio). 


2. 


o-.^0 

«  '■ 
AE 


and 


( a^  ^  822  "  2®21 


AE 


«  1 


T.  «  1 

L®35  "  ®11®'22  ■*■  ®12®2lJ  L  53  AE  J 


and 


IT  a,^  +  agg  -  1  [ 

1^53  "  ®1 1  ®22  +  2®2lJ  L  ^  ^ 


+  ag2  -  (ann22  "  ^12^21^  M 


«  1 


The  method  described  above  is  directly  applicable  to  the  case  where  the 
equations  of  motion  include  3  oigrees  of  freedom.  If  4  degrees  of  freedom  are 
included,  a  similar  technique  can  be  used,  but  must  be  applied  twice — once  to 
get  approximate  factors  for  the  upper  left  3^3  part  of  the  4x4  de  termj.nant, 
and  once  again  to  correct  these  factors. 


B.  APPROXIMATE  FACTORS  FOR  HIGHLY  COUPLED  SYSTEMS 


The  fundamentals  of  this  method  can  be  summarized  briefly  r.s  follows  •  First , 
the  exact  factors  (in  numerical  terms)  must  be  known  for  a  case  vrtiere  the  param¬ 
eters  in  the  equations  of  motion  take  on  typical  values.  Then,  approximate  literal 
factors  (in  terms  of  the  polynomial  coefficients)  are  found  by  solving  the  simulta¬ 
neous  equations  which  relate  factors  and  polynomial  coefficients;  in  this  process, 
numerically  small  terms  hc  -v  rici'X  be  .nool..  ted.  5.u<  n,  because  the  polynomial 
coefficients  are  defined  in  terms  of  the  stability  derivaWi  veo  in  the  equt\:,ions 
of  motion,  the  approximate  factors  can  also  be  expressed  in  these  terms.  A  more 
detailed  description  of  this  method  is  best  presented  in  the  form  of  a  set  of 
instructions.  Although  a  transfer  function  denominator  is  considered  in  the 
following  set  of  instz’uctions,  the  method  is  also  applicable  for  finding  numerator 
factor!' ; 

I .  obtain  ej:act  transfer  function  factors  in  ni;merical  terms  for  a  typical 
set  of  parameter  values 
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2.  assume  the  denaninator  of  the  transfer  function  factors  as  (e.s«;  Tor  a 
sixLh-ox-der  denominator) 

8^  +  Es^  +  Cs^  +  rjs^  +  Es^  +  Fs  +  G  =  [s^  +  (2^0))  ^.^s  +  a|p] 

X  [s^  +  (25aj)^^s  +  u|g] 

X  [s^  +  (25<i>)2eS  +  a|J 

3*  e:q)and  the  factors  in  Step  2,  taid  match  coefficients  of  s,  giving  six 
equations  in  six  unknowns  to  solve 

4.  throw  away  those  terms  in  Step  5  that  are  very  small  (by  knowing  exact 
numbers) 

p 

5.  solve  s.'implified  equations  (from  Step  4)  for  (25o)gp)  Ogp,  etc.,  in  terms 
of  B,  C,  D,  E,  P,  and  G.  This  gives  approximate  frequency  and  damping 
terms  as  functions  of  polynomial  coefficients 

6.  expand  the  determinant  of  coefficients  in  the  equations  of  motion  in 
literal  terms 

T.  match  coefficients  of  s  from  Step  6  with  those  in  the  polynomial  on  the 
left  side  of  the  equation  in  Step  2 

8.  throw  away  terms  that  are  very  small  (by  knowing  exact  numbers)  in  Step  ?• 
This  gives  approximate  expressions  for  the  polynoMai  coefficients  in 
terms  of  coefficients  in  the  equations  of  motion 

9-  combine  results  of  Steps  5  and  8  to  get  approximate  expressions  for  fre¬ 
quency  and  damping  terms  as  functions  of  coefficients  in  the  equations  of 
motion. 

The  applicability  of  •  ;.i.<  n  j.c  '"I'ingent  ./n  one  being  able  to  solve  the  cqua 
tions  in  Step  5j  these  equations  proved  to  be  solvable  ir  all  cases  cf'r.sidered  in 
this  ijtuuy. 


APPiiNDIX  C 


DESCRIPTION  OF  CONFIGURATIONS 

The  data  presented  in  this  appendix  describe  the  vehicles  and  the  flight 
conditions  which  were  investigated.  No  data  are  shown  for  Configuration  1 
because  the  study  of  that  vehicle  was  discontinued  early  in  the  project. 

The  information  necessary  for  the  calculation  of  the  aerodynamic  matrices 
may  be  found  in  Table  C-I  (the  aerodynamic  parameters),  and  Table  C-VIII  (the 
flight  conditions) .  Tables  C-II,  C-IIT,  and  C-IV  present  the  mode  shapes  and 
slopes  for  the  three  configurations,  and  Tables  C-V,  C-VI,  and  C-VII  present  the 
C^(f+oo)3  niatricos.  A  profile  view  of  each  airframe  with  the  control  surfaces 
shown  may  be  found  in  Fig.  C-4  through  C-8. 
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TADLB  C-Il 


MODIi  AMD  FHBqfJHHCIEO 

COIlFIOURATiai  2 


MODR  NO. 

) 

2 

4 

FREQUENCY,  '.-pa 

0 

0 

)  J»21 

5-42 

Mole*  Daflcctlorio 

Cn  nn  r.l 

1  .0 

-C  0241 

0.046*i 

Winy 

1 .0 

v/[ 

0.G04i8 

-0.05&3 

Moic  Olopeio 

C'inarrt 

0 

^  .0 

0.000452 

o.oooeS2 

Wing 

0 

1  .0 

-O.OKI'^ 

0.0001 03 

MODB  OllATEO  AlH)  FREQWKUCIiiS 
CONFICUHATIOIl  5 


MODE  NO. 

1 

2 

5 

n  ! 

FNf:QUBNCY,  opo 

0 

0 

mm 

Dofloctiono 

T>ill 

1  .0 

'Jl'j 

-0.0i/>2 

0.150 

1 /h  Chord  Strip  I 

1  .0 

-120 

-0. 00-^76 

-0.0212 

1/4  Chord  Strip  II 

1  .0 

44.  J 

'  .0455 

-0.074'/ 

1/4  Chord  Strip  III 

1  .0 

203 

0.178 

0.0285 

Mode  Siopeo 

Tail 

0 

1  .0 

-0.0865x1 0'5 

0.707x10*^ 

Stream  Slope 

Strip  I 

0 

1  .0 

0.0157x10'^ 

-0.204xlo"^ 

Stream  Slope 

Strip  II 

0 

1  .0 

-3 

0.1»XiO 

0.l68xl0‘^ 

Stream  Slope 

Strip  III 

0 

1  .0 

0.252x10“^ 

0.580x10"^ 

TABUS  C-IV 

MODE  GttAFES  AMD  FREQUENCIES 
CONFIOURATION  H 


MODE  NO.  '  3 


PREQUENCY,  cpo 

0 

0 

o.y48 

1.251 

Mode  Deflectionu 

1/U  Chord  Strip  I 

1 .0 

-50.0 

0.0122 

-0.0256 

1/4  Chord  Strip  II 

1 .0 

25.0 

-0.00585 

-0.0145 

1/4  Chord  Sti'ip  III 

1 .0 

100 

-0.0>*10 

0.0644 

5/4  Chord  Strip  I 

1  .0 

150 

0.01 53 

-0.0477 

5/4  Chord  Strip  II 

1  .0 

175 

0.0241 

0.00212 

5/4  Chord  Strip  III 

1  .0 

200 

0.0632 

0.115 

-25-48  -54-51  45-50  I  -66.75  50-2  555-6 


TABLB  C-VIII 


PllOlfP  COHDITIOND 


CONPIOURATIOH 

DYHAMIC 

PRESOURE, 

pDf 

FOnWAJiD 

OPEBD, 

ft/ooc 

MACil  NO. 

AI/PITUDE, 

ft 

;>,0S)0 

2,060 

1.8»» 

0 

2 

12,960 

2.96 

0 

639 

1,005 

0.97 

20,000 

1  ji97 

1,005 

0.90 

0 

858 

2,450 

2.37 

22,500 

h 

1,717 

2,450 

2.?0 

40,000 

4,250 

2,490 

2.50 

20,000 

Sta. 

0 


Sta. 

250 


Sta. 


I 

I 


900 


Note;  Stations  shown  in  inches 

Area  of  canard  (total  for  toth  sides)  =  7000  in.^ 
Area  of  wing  (total  for  both  sides)  =  6a,400  in.^ 


Figure  C-4.  Configuration  2 
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4 


Figure  C-5-  Configuration  5 


Stations  in  inches 


Elastic  Nw 


Strip  I  Chord  Area  per  Side  |  Span  Sta.jpus.  Sta. 


■1 

180  in. 

It.llfxlO*'^  in.^ 

11c, 

*10^.5 

II 

l)fO  in. 

?. 22x10*^  in.^ 

5*15 

567.6 

III 

100  in. 

2.30x10'^  in.^ 

575 

752.6 

Tail 

100  in. 

2.0  xlO^*'  in. 2 

L 

0 

1100.0 

Figure  C-6.  Aerodynamic  Strips  for  Corfiguration  3 


ATPENDIX  1) 


NUMERICAL  EQUATIONS  OF  MOTION, 

APPROXIMATE  FACTORS  AND  EXACT  FACTORS 

This  appendix  presents  the  ajajority  of  the  numerical  date  for  the  report. 
The  numerical  equations  of  motion  of  the  three  vehicles  are  given  in  Fig.  D-1 
D-2,  and  D-3,  and  Include  Lhe  static  aeroelastic  correction  for  the  modes  not 
included.  Tables  D-I,  D-II,  and  D-III  present  the  exact  factors  for  the  equa 
tions  and  the  values  obtained  with  the  approximation  formulas  presented  in 
Section  III. 
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TABLE  D-I 

NUMERICAL  VALUES  lOR  EXACT  AND  APPROXIMATE  TRANSFER  FUNCTION  FACTORS 

CONFIGURATION  2 


Transfer 

q  =  35.0  psi 

q  =  90.0  psi 

Function 

3  Modes 

4  Modes 

3  Modes 

4  Modes 

Factor 

Approx. 

Exact 

Approx. 

Exact 

Approx. 

Exact 

Appi-ox. 

Exact 

7.31 

7.33 

7.36 

7.40 

22.0 

22.1 

21 .9 

22.4 

(2e<n)sp 

.335; 

.338 

.328 

.336 

.379 

.373 

.358 

.309 

73.5 

73.5 

73.5 

73.5 

61 .7 

61.6 

62.1 

61 .5 

(2Nle 

.074 

.072 

.074 

.080 

.200 

.196 

.204 

.276 

4e 

- 

- 

444 

444 

- 

4l6 

4l4 

(25a))2e 

- 

- 

.305 

.288 

- 

- 

.835 

.790 

l/^Wi 

^  37 

-137 

-125 

4  25 

-264 

-264 

•202 

-201 

4, 

83.0 

83.0 

83.0 

83.0 

88. 7 

88.2 

88.9  : 

;(2^a))„ 

.219 

.219 

.219 

.219 

.299 

.299 

.302: 

.298 

2 

<^2 

- 

- 

452 

451 

- 

- 

457 

437 

(25‘0)w2  : 

- 

- 

.552 

.531  . 

- 

1 .01  . 

1.03 

i/Te, 

.236 

.242. 

.240 

.241 

■.377 

.381 

.377: 

.381 

-I.  ■ 

83.1 

^  85.2 

83.0 

85.1 

88.7 

=  88.2 

89.0  ' 

(25<o)e.  ; 

.08lu 

.084 

.085 

.085. 

.104 

.104 

.107: 

.■104-- 

2 

«^2 

- 

- 

452 

452 

- 

- 

437 

436 

(2?:<o)02 

- 

- 

.207 

.202; 

- 

r 

.482 

.466- 

“fjl 

5-77  ; 

5.85 

5.84  ; 

5.85 

16.2 

16.2 

16.5  - 

16.3 

j(2Ca))|  ,  : 

.370 

.382 

•  373 

.386 

•533 

•551 

.517 

.558 

alp  - 

444 

444 

4l6 

4l6 

532 

(2503)1^2  : 

- 

- 

.086 

.084 

- 

- 

.405 

.381 

'All’  ■ 

- 

- 

-2.43 

-2.39 

- 

- 

-3.99 

-3.62 

- 

- 

2.65 

2.71 

-  - 

- 

4.26 

4.ly 

‘^42 

- 

- 

87.6 

86.5 

- 

- 

101 

95.9 

.0570 


TABLE  D-II 

NUMERICAL  VALUES  FOR  EXACT  AND  APPROXIMATE  TRANSFER  FUNCTION  FACTORS 

CONFIGURATION  3 
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